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Abstract 

Applying the time-dependent variational principle of Balian and Veneroni, we 
derive variational approximations for multi-time correlation functions in field 
theory. We assume first that the initial state is given and characterized by a density 
operator equal to a Gaussian density matrix. Then, we study the more realistic 
situation where only a few expectation values are given at the initial time and we 
perform an optimization with respect to the initial state. We calculate explicitly 
the two-time correlation functions with two and four field operators at equilibrium 
in the symmetric phase. 
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1 Introduction 



The description of the time evolution of a system with initial conditions which are out 
of equilibrium is one of the main challenges in quantum field theory at finite tempera- 
ture. This problem appears in different contexts in cosmology or in particle physics. In 
inflationary models for the primordial universe, the problem is to describe the evolution 
of the infiaton scalar field in order to be abble to predict the duration of the inflationary 
phase. Another important problem is the prediction of the density fluctuations and its 
relation to the distribution of galaxies In particle physics, the experiments at the 
Brookhaven Relativistic Heavy Ion Collider (RHIC) or soon at the CERN Large Hadron 
Collider (LHC) may show a phase transition between a state where the quarks and the 
gluons are deconfined and chiral symmetry restored and the state of ordinary matter. 
The state of the matter formed after the collision is out of equilibrium and the study of 
the dynamics near the phase transition requires a nonperturbative method. The chiral 
phase transition can be studied in a model of scalar fields with an 0(N) symmetry. The 
relaxation of inhomogeneous pions configurations formed during the collision (called dis- 
oriented chiral condensates or DCC) could be the signal of the chiral phase transition 
through a pions distribution with strong correlations 0. In the Standard model, one 
is interested in the dynamical properties of the electroweak phase transition in order to 
describe the baryogenesis. The violation of the baryonic number, which is related to the 
desintegration of non-trivial topological field configurations, the sphalerons, is evaluated 
from a correlation function calculated at high temperature in a classical approximation 

i- 

The scale for the energy density involved in the primordial universe or the length 
scale which is important for the collective modes in Quantum Chromodynamics require 
non-perturbative approaches. Variational methods are therefore particularly adapted. 

In the mean-field Hartree-Fock (HF) approximation for fermions or the Hartree- 
Bogoliubov (HE) approximation for bosons, one minimizes the free energy or, at zero 
temperature, the energy with a Gaussian variational ansatz for the density matrix or the 
pure state. For a time-dependent problem, the corresponding approximations (TDHF or 
TDHB) are obtained by minimizing an action also with a Gaussian variational ansatz. (At 
zero temperature, a time-dependent variational principle was already introduced by Dirac 
0]. It was used by the authors of reference in the context of quantum field theories .) 
These methods are well adapted to the calculation of expectation values. However, they 
are not suitable for multi-time correlation functions. In order to obtain approximations for 
these multi-time correlation functions, a variational principle was recently introduced in 
the context of non-relativistic fermions . The action to minimize has for stationary value 
the generating functional for the correlation functions. This functional is a functional of 
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the time-dependent sources coupled to the operators of interest. The action involves 
two types of variational objects, one A{t) related to the observables of interest and the 
other T>{t) akin to a density matrix. The temporal evolution of A{t) is described by a 
backward equation while that of involves a time running forward. Let us notice 
that this doubling of the number of degrees of freedom with two temporal arrows appears 
already in the formalism of Schwinger and Keldysh for time-dependent problems at finite 
temperature 0. (For a review of the different formalisms used at finite temperature, 
see reference [|]). By restricting the variational spaces for A(t) and V^t) to Gaussian 
operators, one obtains evolution equations which are coupled. The expansion at first 
order in powers of the sources of these equations (or at the second order for the generating 
functional ) provides approximations for the two-time correlation functions. In spite of 
the Gaussian character of the variational ansatze, the approximations obtained in this 
way go beyond the usual mean-field approximations and include correlations between 
particles. This comes from the fact that the variational principle of Balian and Veneroni 
IP] has been precisely constructed to provide directly the quantity of interest, namely the 
generating functional. 

Our aim is to generalize this formalism to quantum field theories. For simplicity, 
we will consider the case of a scalar field described by a $^ interaction in a Minkowski 
space. We will choose sources coupled to local operators such as $(x) or 11 (x) (written 
here in the Schrodinger representation) and to composite operators such as $(x)$(y) or 
$(x)n(y). The second derivative of the generating functional with respect to the sources 
will thus give two-time correlation functions which have up to four field operators. 

In the first part of the paper, we will consider that the initial state is given and char- 
acterized by a density operator -D(to) equal to a Gaussian density matrix. We will write 
the action functional to be minimized (eqs. ( |2.14| ) and ( |2.15D ). Then we will restrict our 
variational spaces to operators which are exponentials of quadratic and linear forms in the 
field operators $(x) and n(x) (called Gaussian operators). An adequate parametrization 
of the variational objects allows us to write the coupled dynamical equations in a compact 
form. We stress again the fact that the variational method we use introduces a doubling 
of the degrees of freedom. We obtain four coupled equations (two for the expectation 
values of one- field operator and two for the expectation values of two-field operators). 
Therefore we have twice as many dynamical equations as the usual mean-field approach. 
The expansion of these equations to lowest zero order in powers of the sources leads to a 
decoupling of these equations whose a half is identical to the time- dependent mean-field 
equations for the $^ theory written in reference We have showed in a previous pa- 
per that these dynamical equations can be put into a classical noncanonical Hamiltonian 
framework with a Poisson structure which is a generalization of the standard Poisson 
bracket [IC]. At the lowest order we thus obtain an approximation for the single-time 
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expectation values of one or two field operators. The expansion of the coupled equations 
to first order in powers of the sources provides an approximation for the two-time cor- 
relation functions with two, three or four field operators. Their evaluation proceeds in 
the following way : first one solves the mean-field TDHB equations, then one calculates 
the kernel which corresponds to small deviations arround the mean-field solution (which 
is analogous to the random phase approximation (RPA) kernel of reference |^) and fi- 
nally one solves a backward dynamical equation which involves the kernel calculated for 
the mean-field solution. These results have been obtained for the generating functional 
of correlation functions with time-ordered product of the field operators. We also give 
the corresponding expressions for the retarded Green functions. Indeed, these retarded 
Green functions appear when one studies the response of the system to a small external 
perturbation. A less complete version of this first part of the paper has been published 



in reference 11 



In section 2, we have supposed that the initial state -D(to) was known. In practice, 
this is not the case. Only the expectation values of a few observables are known. In section 
3, we use the form of the variational principle introduced by Balian and Veneroni which 
optimizes also the initial state 0. The approximations obtained in this way for the two- 
time correlation functions include not only the correlations between particles originating 
from the dynamics but also those present in the initial state. One introduces again 
two variational objects which now evolve in imaginary time. We study the symmetric 
phase. We will consider the two-time function with two field operators $ and the two- 
time function with four field operators $. We will relate it to the standard polarization 
function 11;^. At equilibrium, in the symmetric phase, the case of the correlation function 
with two field operators is special since the approximations obtained with or without 
optimization with respect to the initial state coincide. For the correlation function with 
four field operators, we will be obliged to combine the variational approximation with a 
perturbative approximation valid for small values of the coupling constant, by neglecting 
or taking into account at first order a loop term in the dynamical equations. However, the 
solutions remain nonperturbative since they depend on the self-consistent HB solution. 
At the lowest order, by neglecting the loop term, the two-time function with four field 
operators is simply related to the two-time function with two field operators according 
to the Wick theorem. At equilibrium, in the symmetric phase, we give to first order 
the expression for the two-time polarization functions with four field operators : the 
retarded Green function n/j(q^, t', t") (eq. 4.133a) and the non-retarded Green function 
n(q^,t',t") (eq. ( p.l44| )). We verify explicitly that these functions depend only on the 
time difference t' — t" and not on the initial time to. This is not the case for the non- 
retarded Green function if, like in the section 2, we optimize only the dynamics and not the 
initial state. We check also that our approximations for the two-time correlation functions 
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satisfy the fluctuation-dissipation theorem. The last section contains some conclusions and 
perspectives for future work. 

2 Variational approximation for the two-time corre- 
lation functions : optimization of the dynamics 

2.1 Definition of the generating functional for correlation func- 
tions 

We define the operator A{J,t) which is a functional of the sources J and a function of 
the time t in the following way : 

A{J,t) ^T&L^i dt' y d^x J^{^,t')^^{x,t',t) + J^{x,t')U^{x,t',t) 

+ 1 d'x d'y J*n(x, y, t') ($^(x, t', t)U^{y, t\ t) + n^(y, t' , i)$^(x, t' , t)) , 

+ J**(x, y, t') $^(x, t', i)$^(y, t', t) + J™(x, y, t') n^(x, t' , t)U"{y, t', t)} 

(2.1) 

where T is the time-ordered product of operators and $^(x, i',i) and U.^{yi,t',t) are 
respectively the field and momentum operators defined in the Heisenberg representation. 
They satisfy the boundary condition : $^(x, t) = $(x) and U."{x,t,t) = n(x), $(x) 
and n(x) being the operators in the Schrodinger representation. We use the short hand 
notation 

A{J, t) = Texp (^i dt' Jj{t') Oj{t', t)^ (2.2) 

The sources Jj{t') (which are local or depend on two space points) are turned on between 
t' — t and t' — +00. 

The functional A{J,t) satisfies the equation 

j^A{J, t) = i[AiJ, t),H]-i AiJ, t) Y: Jj{t)Oj (2.3) 

with the boundary condition A(t = +oo) = 1. We will consider the case of a self- 
interacting scalar field in a Minkowski metric described by the Hamiltonian H = J d^x7i{x) 
with : 

1 1 -» Tn^ h 

H(x) = -n2(x) + -[V<l>(x)]^ + ^$2(x) + -$^(x) . (2.4) 

We work in three spatial dimensions. The constants mo and b are respectively the bare 
mass and the bare coupling constant. 

The generating functional for the causal Green functions writes 

Z{J,to)=Tr{D{to)A{J,to)) , (2.5) 
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D{to) being the initial state. The statistical operator D(to) may represent a thermal 
equilibrium or a nonequilibrium state. At zero temperature, it reduces to a projection 
operator. In addition to the sources J, the generating functional Z depends on the initial 
time to. This is a difference with the generating functional usually considered in quantum 



field theories 12 . We want to evaluate the functional derivative of 



W{J,to) = -t \nZ{J,to] 



(2.6) 



W{J,to) is the generating functional for the connected Green functions. Let us note that 
for hermitic operators A{t) and D{t), W{J,t) is real. The expansion in powers of the 
sources of W{J,t) writes : 

W{J, to) = -i In(no) + dt' { J d^x J$(x, t') y?(x, t') + J d^x d^y J$$(x, y, t') G(x, y, t') + ...} 

+ ^ y ^ " dt' dt" { J d^Xi d''x2 M^l, f) J$(X2, t") C|^(X1, X2, t\ t") + ...} 



+... 



(2.7) 

where uq = TrD{to) ( we choose not to normalize the operator D{to)). We introduce the 
following notations for the expectation values of one and two field operators : 



^(x,t) = — Tr ($^(x,t,to)I^(tc 
no ^ 

^(x, y, t) = — Tr ($^(x, t, to) $^(y, t, to) D{to)) - v^(x, t) ^{y, t) . 
no ^ ^ 

The two-time causal functions with two field operators are defined according to 
CU^, y, t', t") = ^Jr{T <|.^(x, t', to) $^(y, t" , to) /^(to)) - ^(x, t') ^{y, t") . 
Similarly, we define the functions C|n, C^n- From the expansion (|2.7| ), we have 



(2.8) 
(2.9) 

(2.10) 



6^W 



i 5Jq>{x,t')5J^{y,t"] 



j=o 



C|.$(x, y,t ,t 



(2.11) 



We define also the three-point and four point two-time correlation function et : 
6^W 



C\^,y,z,t',f 



5J$(x, f)5J$$(y, z, f ^ 
+^(y , f ) (C2(x, z, t', f ) + Vi^, t') ip{z, t"))+ip{z, r ) {C^ (x, y , t' , t" ) + y.(x, t') ip{y, f 

(2.12) 

T jn ^ ^C\^,y,z,n,t',t") 

()J$$(x,y,t')()J$$(z,u,t ) j=o 

+ (C'i^, z, t', t" ) + ^(x, t') ifiz, t")) {C'iy, u, t', t") + ifiy, t') ^(u, t" ^^.13) 
+ (C2(x, u, t', r ) + <^(x, t') <^(u, t" )) (C2(y , z, t', t" ) + ^(y , t') <^(z, f ^ 
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2.2 Variational evaluation of the generating functional 



We will use the time-dependent variational principle of Balian and Veneroni [^] to obtain 
an approximation for our quantity of interest Tr{D{to)A{J,to)). This principle is an 
illustration of the general method developped in reference [l^ to construct a variational 
principle suited to the evaluation of a desired quantity. The basic idea is to consider 
the equation ( p. 3D which defines A(J, tg) as a constraint and to introduce a lagrangian 
multiplier T'(to) associated to this constraint. 

We define the action functional : 

Z {A{t),V{t)) = Tr {A{to) V{to)) + Z,y^ , (2.14) 

where 



Zdyn = Tr dt V{t) - I \A{t), I d'xHi^)^ + tA{t) Mt)Qj 



(2.15) 



We have written in a compact form the term which involves the sources J : 



Ei Jj{t)Qj = J (j*(x, t)<i>(x) + jn(x, t)n(x) 

+ J d\ d'y J*"(x, y, t) (<l'(x)n(y) + n(y)<l>(x)) (2.16) 
+ J d'xd'y (J**(x,y,t)$(x)<l.(y) + J™(x^y^^)n(x)n(y); 



The variational quantities of the action functional Z are the observable-like operator 
A(t) and the density-like operator T>{t). We look for the stationarity of Z with respect to 
variations of A(t) and Vit) with the boundary conditions : A(t = +oo) = A(t = +oo) = 1 
and 'D{to) = D{to)- D{to) is the initial statistical operator which we assume to be given 
and equal to a Gaussian density matrix. By construction of the variational principle, 
the generating functional for the connected Green functions will be approximated by the 
value of the functional Z at the stationary point : 

l^(J,to) = -nnZ,i . (2.17) 



2.3 Choice of the variational spaces 

We restrict ourselves to trial operators A{t) and 'D{t) which are exponentials of quadratic 
and linear forms of the field operators $(x) and n(x) which we shall loosely call Gaussian 
operators. 

A Gaussian state 'D{t) is completely characterized by the vector a(x, t) and the 
matrix S(x, y,t) defined by : 

«(x,*)=f ) . (2^18) 

V-Z7r(x,t) ) 
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/2G(x,y,t) -zT(x,y,t)\ 

^(x,y,t) = , (2.19) 
V-«T(y,x,t) -25(x,y,t)/ 

where 

(^(x,t) = -^Tr(P(t)$(x)) , (2.20) 

7r(x,t) = ^Tr(P(t)n(x)) , (2.21) 

n{t) 

G(x, y, t) = 4t 2^K^^W3'(x)$(y)) , (2.22) 
n{t) 

5(x, y, t) = ^ Tr(P(t)n(x)n(y)) , (2.23) 
n{t) 

T(x,y,t) = -^Tr(P(t)(<l(x)n(y)+n(y)<l(x))) , (2.24) 

n{t) = TrVit) , (2.25) 

and, for any operator O : = — Tr{T>{t)0). The hermiticity of 'D{t) imphes that the 
matrices G, S and T are real and that G and S are symmetric in the space indices. 

The Heisenberg invariant is equal to : 
/(x, y) = / d^z (4 < <l(x)<l(z) X n(z)n(y) > 

- < $(x)n(z) + n(z)$(x) > < $(z)n(y) + n(y)$(z) >) ^ ' ' 

/(x,y) = j (-Hn(x,z)S22(z,y) + Hi2(x,z)Si2(z,y)) . (2.27) 
For a pure state : 

Sn S22 - + 1 = (2.28) 

/(x,y) = 53(x-y) (2.29) 

In Appendix A, we give the relations between our parametrization of a Gaussian state 
and the parametrization used by the authors of reference 



We then introduce the operators % = A{t) V{t) and % = V{t) A{t) 0. The 

operator 7^ is a Gaussian operator which is characterized by the following expectation 
values : 

n{t) = Tr% = Tr% , (2.30) 

^,(x,t) = -lyTr(T,<l.(x)) , (2.31) 

7rfe(x,t) = ^Tr(rfcn(x)) , (2.32) 

G,(x, y^i) = ^ Tr(T,<l(x)<l(y)) , (2.33) 
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^,(x, y,t) = ^ Tr(T,n(x)n(y)) , (2.34) 

r,(x, y^i) = ^ Tr(rb($(x)n(y) + n(y)<|.(x))) , (2.35) 

with i>(x) = $(x) - (pbi^,t) and n(x) = n(x) - 7r;,(x, t). 

Similarly the operator T^is a Gaussian operator which is characterized by the quan- 
tities n, (fc, TTc, Gc, Sc, Tf, ( obtained by replacing 7^ by 7^ in (p.31|) - (|2.35|) ). For Vit) = 1, 



one has ipb{'x,t) = {pc{'x,t) = (pa{'x,t) (and similar definitions for TTa,Ga, Sa,Ta). For 
A(t) = 1, one has (pb{^, t) = V5c(x, t) = V5d(x, t) (and similar definitions for tt^, Gd, Sd, Td). 
In Appendix C, we give the expressions of 0:5, Sf, and Sc as functions of the independent 
variational quantities eta, and a^, which characterize respectively A and V. 

With this choice of the trial spaces for A{t) and T>{t) , the Wick theorem allows us 
to express the functional Z in the form : 



Z= n(to) + / dt —\vit)=cte - i n{t) j d^x {£c{^,t) - £b{y^,t)) 

Jtu dt J ^ (2.36) 



+in{t) J d^Xid^X2 JCc{xi,X2,t) 



where the energy density £^(x, t) = Tr {T>{t) H(x, t)) is given, for the Hamiltonian density 
D,by : 

£(x, t)= \ vr2(x, t) + \ (V¥.(x, t)f + ^ ^^(x, t) + A y,^(x, t) 

+ ^ 5(x, X, t) - i AG(x, y, t) + ^ G(x, x, t) (2.37) 
+ -G2(x,x,t) + ^y.2(^^^)^(^^^^^) 



and ^ 

In n{t) = lnna + In Ud - -at-cTaa-d 



+ ^ Indet (ai(l + r)(H„ + Hd)-^(1 + r 



(2.3^ 



2 

where ab-c = o^b — ac and r and cxi are the following 2x2 matrices : 









^ 1\ 
















v-1 0/ 



^1 = L J > ^ = . J • (2-39) 



The last term of the functional ( ^.3(i| ), = Tr [% J2j Qjj = I d^xid^X2lCc{^i, X2, t), 
involves the sources J and the expectation values indexed by c. 

With the parametrization we use, the boundary conditions T>{to) = D{to) and 
A{t = +00) = 1 correspond respectively : 

nd{to) = no , adito) = , Ed{to) = Eq , (2.40) 
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where no, ao and Sq characterize the Gaussian initial state D{to), and 

Uait = +00) = 1 , E-\t = +00) = , ^aa{t = +00) = . (2.41) 

' — 'a 

We will write the conditions for the stationarity of the functional Z with respect to the 
variations of n^, a^, and n^, a^, S^. 



2.4 Dynamical equations with sources 



By varying the expression ( p.36|) with respect to n^, a^, S^, with the boundary conditions 



ridito) = no , adito) = ao , Ed{to) = Hq , (2.42) 



we obtain the evolution equation for Ua, aa and S^. Integrating ( |2.36|) by parts and varying 
with respect to n^, with the boundary conditions 

na{t = +00) = 1 , E-\t = +00) = , ^aa{t = +00) = , (2.43) 

' — 'a 

we obtain the evolution equations for Ud, otd and S^. In general the evolution equations 
for n^, Od, '^d and those for n^, are coupled. The solutions n^, ad, and n^, 

depend on the sources. Owing to this dependence, we will he abble to obtain non trivial 
approximations for the correlation functions in spite of the Gaussian nature of the trial 
operators. 

By combining the evolution equation for nd,ad, with those for n^, eta, Sa, the 
dynamical equations for the expectation values at, "^b and Oc, '^c can be written in the 
following compact form : 

i ctb = T Wb , (2.44) 

iac = T {wc- w^) , (2.45) 

tEb = 2{Ebnbr-rnbEb) , (2.46) 

zE, = 2 (S, (K - X^) r - r {H, - J^) H,) . (2.47) 

We have also n{t) =0. 

The vector and the matrix Tib are defined from the variation oi < H >b= 
TriTbH) : 

6<H>b= J d^x Wb{^, t) 5afe(x, t) - \ j d^x d^y tr [7^b(x, y, t) 6Eb{y, x, t)] . (2.48) 

Similarly, w^'Hc and iu|^,X^ are defined from the variation of < if >c= Tr{TcH) and 
from the variation of the source term Kc- 
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Therefore, by using a convenient parametrization for the two variational objects 
T>{t) and A(t), we have been abble to obtain dynamical equations in a very compact form 
even for finite values of the sources. The dynamical equations for ad, and aa, have 
a more complicated form. They are given in Appendix D where we give also the explicit 
expressions of the vectors w, w^^ and of the matrices Ti and in the case of the $^ 
theory. In spite of their simple form, the equations ( p.44D -( p.47[ ) are coupled because the 



solutions ah, S;, and Oc, 2c do not satisfy simple boundary conditions. 

Let us note that the coupled equations for ad, {aa, S^) do not conserve neither 
the energy < H >d {< H >a) nor the entropy Sd {So) (see Appendix B for the definition 
of the von-Neuman entropy for bosons). On the other hand, the dynamical equations for 
a^ and (eqs. ( p. 44 ) and ( 2.46 ) ) which have the same form as the mean-field TDHB 



equations conserve the pseudo-energy < H >b and the pseudo-entropy Sb- 



The expansion in powers of the sources of the stationarity conditions ( p.44| )- ( |2.47| ) 



will provide approximate dynamical equations for the expectation values and the corre- 
lation functions defined in Eqs. ( p78|) -( pyi3l) . 



2.5 Expansion in powers of the sources 
2.5.1 TDHB equations for the expectation values 

We will use the upper index^^^ for the solutions of the dynamical equations when there 
are no sources. The limit with vanishing sources corresponds to 

«r^ = «f = «r , Hr = Hf=sir (2.49) 
From equations ( |(J.4| ) and ( |(J.5| ) and from 

ab - ac = 2 T ^ ("a - ttd) , (2.50) 

~r ^d 

we obtain that in the limit with vanishing sources 

S-^ = , ^aa = 0. (2.51) 

' — 'a 

We have also Wb'-^^ = Wc^°^ = and Hb^^^ = Uc^^^ = n^°\ In this limit, the dynamical 

equations (p.2| ) and (p.3|) for ad and Ed become (we suppress the index d) : 

^ d(o) = r , (2.52) 

, = - [(SW + r) - r) - - r) + r)] . (2.53) 
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These equations are the analog for the $^ theory of the time-dependent Hartree-Bogohubov 
(TDHB) equations in nonrelativistic physics. They are equivalent to the dynamical equa- 
tions obtained in reference [^] where the authors used an alternative form of the Balian- 
Veneroni variational principle suited to the evaluation of single-time expectation values 



In reference |jlO|, we have shown that these TDHB equations for the $^ theory 
can be written as classical dynamical equations with a nonsymplectic Poisson structure 
with the free-energy density playing the role of a classical Hamiltonian. The Heisenberg 
invariant (p.26|) appears as a structural invariant of the Poisson structure. 

The static HB solution satisfies Wq = that is, in the uniform case : 

b b 

{ml + -ipo + -Go(x, x))(^o = (2.54) 
2 

and the so-called gap 

So = -r cothipHoT) (2.55) 



/3Ho = irlnf=^^') . (2.56) 



which can also be written as : 



2 VSo-r 

More explicitly, in the uniform case, by going to the momentum space : 

/ J- coth ^ \ 



and for the $ theory, tOp = y—go{p) with 

^o(p) = -ip' + ml+ ^-ipl + ^Go(x, x)) . (2.58) 

The solution y^o = of equation (|2.54| ) corresponds to the symmetric phase while 
the solution y^o 7^ corresponds to the broken phase. 

We define the self-consistent mass m{(3) depending on the temperature as the solu- 
tion of the gap equation : 

m\l3) =ml + + ^Go(x, x) (2.59) 
By using eq. ( ^.54| ), we have also : 

m'{P)=c(ml + ^Go{^,^)^ (2.60) 
with c = 1 in the symmetric phase and c = —2 in the broken phase. 
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The gap equation ( p.59| ) has ultraviolet divergences and requires a regularization. 
The divergent part Gq'^(x, x) of Go(x, x) comes from the zero temperature contribution. 
By introducing a cutoff A in momentum, 



G^fx, x) 



87r2 



- log ( ) 

em J 



(2.61) 



where m = m(T = 0). The gap equation can be put in a finite form by choosing the 
following renormalization conditions [1^ : 



2 I fcA^ 

.2 _ '"o + Te^ 



^/^ = 1 , X > ^ = 0,±1 (2.62) 

1 2 1 , / 2A \ , , 

where /i is an arbitrary mass scale. The renormalized gap equation writes at zero tem- 
perature : 

m' = ef.' + '44^'^og(%] (2.64) 



In the next sections, the dynamical equations for the two-time correlation functions 
will involve the bare coupling constant b. This one will be understood as a function of 
the cutoff A. We will solve the dynamical equations for a finite value of the cutoff. In 
the static case, we will show that the final formula obtained for the two-time correlation 
functions can be written with the renormalized parameters. 

At finite temperature, we can isolate the divergent part of Gq{'x,x) according to : 

Go(x,x) = / - — h / 7 — 2 (2.65) 

Therefore, the gap equation at finite temperature writes : 



2 , 9Rif^) 2/^M fm'^iPy 



m\/3)= e/.^ + ^m^(/5)lo , 
167r^ y /i^ 



(2.66) 



167r2 ''^^ J (27r)3 el^^>= - 1 
2.5.2 Two-time correlation functions : 

The first derivatives of iy(J, to) = —i In^si with respect to the sources are equal to the 
expectations values with the index c. Indeed, the functional Z depend on the sources 
both explicitly and implicitly since the approximate solutions for V{t) and A{t) depend 
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on the sources. However, at the stationary point, only the exphcit dependence contributes 
to the first derivative : 

dZ 5Z 



dJj{t) 5Jj{t) 

which gives for instance : 



tTr{V{t)A{t)Q,) , (2.67) 



The expressions for the second derivatives of W are much more comphcated. The intro- 
duction of sources coupled to the composite operators $(x)$(y), $(x)n(y) and n(x)n(y) 
together with eq. ( 2.68 ) allows us to obtain dynamical equations for the two-time cor- 
relation functions with three or four field operators merely from the expansion of the 
expectation values ac and at the first order in powers of the sources. 

From the first order corrections ac — a^^^ and Sc — H*^*^-*, we define the two-time 
correlation functions (3 and S (/? is a vector and E is a matrix) : 

ae(x,t)- a^'^\^,t)c^i j'^ dt" |y'rf3^/3*(x,y,t,t") J*(y,0+/5"</" 

+ j d'x,d'x2 /3**(x, xi, X2, t, t") J**(xi, X2, t") + + /5™ J™| 

(2.69) 

S,(x, y, t)- SW(x, y, t) ^ ^ dt" | J d'x, S*(x, y, Xi, t, t") J*(xi, t") + 

+ j d'x,d^X2 S^x, y, xi, X2, t, t") J**(xi, X2, t") + S*n + S™ J™} 

(2.70) 

We have : 

/3f *(x, y, z, t, t" ) = ^ (x, y, z, t, t" ) . (2.71) 

The functions (3 and S provide approximations for the exact two-time correlation functions 
C2, and defined by Eqs. (|l0|)-(|ll : 

C|^(x,y,t,t")^/3f(x,y,t,f'), (2.72) 

C3(x,y,z,t,t")^ /5r(x,y,z,t,f' ) 

-^(°)(y,f') (/3f(x,z,t,f') + ^W(x,t)^W(z,f')) , (2.73) 

VHz,t") (/3f(x,y,t,t") + ¥.(°Hx,t)¥^(°Hy,t'^ ' 



C4(x, y, z, u, t, t" ) ~ ^ Sff (x, y, z, u, t, t" ) 

- r/5f(x,z,t,t") + ^^°nx,t) ¥^(°nz,t")) (/5f(y,u,t,t") + ¥'^°ny,^) ^^°nu,t") 

- (/5f (x, u, t, t") + V. W (x, t) ^ W (u, t" )) (/3f (y , z, t, t") + V'^^^ (y , t) ^^'^ (z, t" ) 

(2.74) 

The expansion in powers of the sources of the dynamical equations (|2.45| ) and (|2.47| ) 
for ac and Sc yields : 

i 5ac = T {6wc — Sw'^) , (2.75) 
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i6E^ = 2 



(2.76) 



(SHc - 61],) r-r ((5K - 5Tk) 

In these equations the matrix Ti. has to be evaluated for the TDHB solutions a^^\ S^^^ of 
eqs. (|23^) and ( ^331) . The variations of ti;^(x,t) (eqs. ( pl8D -( pl9D ) and of J^(x,y,t) 
(eqs. ( L).2(J[ )-( p.22D ) with respect to J{t") will give terms proportional to 6{t — t"). From 
the variations 6w'^ and 5?-^^, we define the matrices tij,Tijk,rij,k, Rij,ki (which are the 
analogues of the random phase approximation (RPA) kernel of Balian and Veneroni 0) : 

5<(x, t)= J d?y ti,(x, y, t) (5«^^(y, t) j d^y d^z Tij^i^, y, z, t) SElj{z, y, t) , (2.77) 

y^t) = J "^ijA^^ y> z, t) ^afe(z, ^)- \ j Rij,kii^^ y> z, u, t) (5S;=fc(u, z, t) . 

(2.78) 

These matrices depend on the TDHB solutions a^^^ and 'E^^\ Their expressions in the case 
of the <1>^ theory are given in Appendix D. With these notations, the dynamical equations 
for the two-time and two-point correlation function /3* and the two-time and three-point 
correlation function write : 



^|/3*(x,y,t,f') = r 



t(x,z,t) /3*(z,y,t,t") - - T(jfc(x,z,u,t) Sfj.)(u, z, y, t, t") 
/I' 



+1 







z|E*(x,y,Xi,t,r 



2 

+ 



5(x-y) 5(t-f' 



(2.79) 
(2.80) 



We obtain similar equations for the three-point functions f3'^^, f3^^ and the four-point 
functions S**, and S™. 

These dynamical equations are not sufficient since the boundary conditions are 
H^^(t = +oo) = 0, i^ctait = +oo) = and ad(x, to) = ao(x), Hd(x, y, to) = So(x, y). But, 
at this step, we have used the expansion of only one half of our variational parameters. 
We will use also the dynamical equations for the two-time functions m and I defined from 
the expansion in powers of the sources of S'^aa and of at the first order : 



^aa{t)c^i dt" m\r,t)Jk{t") 



'to 
r+oo 

s-i(t)~z / dri''{t",t)Jk{tn . 

Jto 

(We have ommitted the space indices) 



(2.81) 
(2.82) 
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From the relations of Appendix C, we will be abble to obtain relations between 
the functions T,''{to,t") and l^(t",to) one the one hand and the functions l3^{tQ,t") and 
m'^{t",to) — I'^it" ,to)ao on the other hand. Indeed, that at the first order : 



6E 



-—6E — 



From eqs. (| 



and (|C.2|) , we obtain at the first order : 



(2.83) 



(2.84) 



and 



6ar 



- r) 6{: 



an 



:(o) 



-r ^d 



^ ^ 

' — 'a ~l~ ' — 'd • — 'a ~l~ ' — 'd 



+ Sad ■ (2.85) 



At the initial time, = 0, Sa^ = and we have the following relations between the 
two-time functions : 



E^'(i^o,t") = -(2o-r) I'it^to) (So + r) 
/3\to,f') = (So - r) im'if\to) - /'(t",to)«o) 
By using the dynamical equation for S^^^ at the first order 



(2.86) 
(2.87) 

(2.88) 



and the variations 6wc — Swt and STic — STib, we obtain a backward dynamical equation 
for l^{t",t) : 



t f/{t" , t),, = 2 r^,- TH (m^ - /^a(°)); - i?,,- / r - r Z S 



(0) 



+ 



5Jt 



5(t-f' ) 



Zfc 



(2.89) 



where the matrices r and R have been defined in (|2.78| ). 

Similarly, from the dynamical equation for S^^Oa at the first order, we obtain a 
backward dynamical equation for m^{t" , t). It is actually more convenient to consider the 
function m^{t" ,t) — l^{t" ,t)a^^\t) which is solution of : 



i| (m*^(t",t) -/'^(t" ,t)aW(t)) = tT {m\t" ,t) - l\t\t)a^'^\t) 

1 
2 



2 oJk 



(2.90) 

where the matrices t and T have been defined in ( p.77| ). The reason to use this function 
will appear with the introduction of the function F{t',t") bellow (see |2.100|) . 
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The boundary conditions for these two equations are 



= 0, m(t",t) = pour t > t" (2.91) 

/^(t",^) = , ) = (2.92) 

/ff (x, y, z, u, t" , t" ) = -^(5^(x - z)53(y - u) + ^^(x - u)5^(y - z)) (2.93) 

lf^{^,tn = 2^ , C(f',f') = 2^ , l^^{t\r) = l (2.94) 
The other elements of the matrix I are equal to zero. 

mf(x,y,t",t") = -^5'(x-y) , mn(x, y, t", t") = - V(x - y) (2.95) 

mr(x,xi,X2,t'',r)-Zf*(t",t")af (t'') = -<^^°Hx2,t")5'(x-xi) + sym. (2.96) 

mfn(x,xi,X2,t",t")-/fn(t",t")af^(t") = -^^°Hx2,t")5'(x-xi) + sym. (2.97) 

m|n(x, xi, X2, t" , t" ) - Itl'it" , )«r ) = -^""^ (X2, )'^'(x - xi) + sym (2.98) 

m™(x,xi,X2,t",t")-C(f' ,f' )«f^(^") = -^7r(°)(x2,t")(53(x-xi) + sym (2.99) 

where sym. means the same terms obtained by exchanging X2,x and xi. 

The same matrices t, T, r, R coming from the expansion of < > at the second 
order arround the TDHB solution appear in equations (|2.79|) ,( p780| ) for /? and S and in 
equations (|2.89|) , (|2.9CI| ) for m and /. This leads to consider the following quantity : 

F\t',t") = 2 {m\t',t) - a^''\t)l\t',t)) (3\t,r) -^tr {l\t' ,t)J:\t,t")) (2.100) 

where m is the transpozed vector of m. More explicitly, by reintroducing the space 
coordinates, F'^{t',t") writes in the case of the functions m'^ , l'^ , jS"^ and : 

F*(t',t")= 2 J £x[m^((K,^i,t',t)- j c/3Ma(°)(u,t)/*(u,x,Xi,t',t)^/?*(xi,X2,t,t") 
d^x d^ztr (/*(x,z,xi,t',t)S*(z,x, X2,t,t")) 

(2.101) 

The quantity F'' is independent of the time t : from the dynamical equations ( |2.79| ),( p.80| ) 
and ( |2.89|) ,( pr90D , one shows that its derivative with respect to t is equal to a sum of 
functions 6{t — t') and 6{t — t") : 



^fF'' = % ^\t,f)^6{t - t') + i2{m\t',t) - &^'^\t)l\t' ,t))T^6{t - t") 
oJk oJk 

-\tr {^^S(t,t")5(t - t') + U\t\m^'\t)^-ff-r5{t - t") - tl\t',t)T^-^E^'\mt - t") 

(2.102) 
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The fact that does not depend on the time t will save us to solve one dynamical 
equation. 

By considering the local sources J*(x, t) and by integrating equation ( p^.l02| ) from 
t = +00 to t = to, we obtain an approximation for the two-time and two-point function. 
For t" > t' : 

f3f{t',t")= 2 



1 



^2.103) 



(2.104) 



/*(t',to)S*(to,t")] 
In the uniform case, we have in momentum space : 

Pf{p,t',t") = 2(m*(-p,t',to)-«o/*(0,p,-p,t',to))/?*(p,to,0 

-- J d'k tr [/*(p + k, -k, -p, t', to)S*(k, -p - k, p, to, t") 

Our conventions for the Fourrier transforms are defined in Appendix D. 

These two last expressions will still be valid when we will optimize with respect to 
the initial state in section 3. On the contrary, the following expression will not be valid 
in section 3. By using ( p.86|) and (|2.871) , we obtain : 

/?f(t',t") = 2(m*(t',to) -«o/*(t',to)) (Ho-r) (m*(t", to) - to)«o) 
+Ur [f{t',to){^o-r)f{r,to){^o + r) 

For t' > t" , one interchanges the operators (f)^{t',to) and (j)^{t" ,to). By using the expres- 
sion of the retarded Green functions (see bellow), we show : 

Pf{t',t") = 2(m*(t",to) -ao/*(t",to)) (5o - r) (m*(t', to) - /*(t', to)«o) 
"/*(t",to)(Ho-r)/*(t',to)(Ho + r)" 



(2.105) 



1 

+-tr 
2 



. (2.106) 



For the correlation functions with the ant i-T-pro duct, it is sufficient to replace in 
the functional Z ( |2.36| ) the source term JCc by the analogous term /C^. This is equivalent 
to exchange the matrices Sq — r and Sq -|- r in the previous formula (or to exchange t' 
and t"). 

From the function F'^(t',t"), we obtain also, by considering now the bilocal sources 
J**(x, y, t), an expression which relates the approximations for the two-time and three- 
point function /5f*(t',t") and for the two-time and four-point function Sf*(t',t"). For 
t" > t' : 

/3r(xi, x;, x'2, t', t")^(°)(x2, t') + /5r(x2, xl, x'2, t', t")v5(°)(xi, t') 

(Sff (xi, X2, x;, x'2, t', t") + Eff (X2, xi, x;, x'2, t', t")) = 
-2 Jd'x (m**(x,Xi,X2,t',to)- y'rf3^ao(u)/**(u,x,Xi,X2,t',to))/5**(x,x;,x'2,to,t"; 
+ ^ y" d^xd^ztr /**(x, z,Xi,X2,t',to)S**(z,x, x^,X2,to,t") 

(2.107) 
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In the uniform case, we have in momentum space : 

'PfiPi, q, -(Pi + q), t', t")6%p2) + /5r(p2, q, -(P2 + q), t', t")S%p,)] ^^'\t') 
Sn*(pi, P2, q, -(Pi + P2 + q), i', t") + Sff (p2, Pi, q, -(pi + P2 + q), t', t"] 



Pi - P2, Pi, P2, t', to) - ao^**(0, -Pl - P2, Pi, P2, t', to] 



4 
-2 

x/?**(pi + p2, q, -q - Pl - P2, to, t") 
1 



d-^ktr 



I (k - Pl - p2, -k,pi,p2,t',to)S (k, -k + pi + p2,q, -q-pi -p2,to,t") 

(2.108) 



Let us recall that the function is related to /?f*(t',t") by : /?f* = . By 



comparing the dynamical equations for /5f*(t',t") and S* (t',t") we obtain on the other 



hand : — ^12 + ^21 expression for E22 in terms of Pf and two integrales 

over the momenta which involve Sf* and = 2/5f*. The functions are there- 
fore completly known if /5f* and Sf* are known. However, we note a difficulty in the 
asymmetric phase since we have two equations ( |2.103|) and (p.l07| ) for the three functions 
(3f{t',r),(3f'^{t',t") and Sff (t',t"). 

In the symmetric phase, equations ( p.l05| ) and ( p.l07| ) simplify. Indeed, from 
the boundary condition l^{t\t') = and the dynamical equation ( |2.89| ) for /, we deduce 
/*(f,t) = for t < t'. Similarly, m**(t',t) - = for t < t'. Therefore, we 

have in the symmetric phase and for t" > t' : 

/3f (p, t', t") = 2 m*(p, t\ to) (Ho(p) - r) m*(p, t\ to) (2.109) 

and for the four-point function : 



Sfi (Pl, P2, P3, - E Pi, t', t" 

By using ( p.86| ) : 

S?f(pi,P2,P3,-Ep.,t',t") 



/ 



(27r)3 



tr 



I (P-Pl -P2,-P,Pl,P2,t',to) 



tr 



Sf*(P, -P + Pl + P2, P3, - Yl Pi' ^0, t") 

(2.110) 



/ (P - Pl - P2, -p,Pl,P2,t',to)(So(p) - t) 



^**(P, -P + Pl +P2,P3,-I]Pi,t",to)(2o(-p + Pl +P2) +7-) 

(2.111) 

Let us emphasize that it is Hg, the static HB solution and not S*^*^) which appears in 
expression ( |2.111D . 

For t' > t" , we have, still in the symmetric phase : 

Pf{p,t',t") = 2m*(p,t',to) (So(p) + r) m*(p,t",to) (2.112) 



s?f(pi,P2,P3,-Ep.,t',n 



-tr 



I (P - Pl - P2, -p,pi,P2,t",to)(2o(p) - r) 



^**(P, -p + Pl + P2,P3, -X1p»'^'''^o)(2o(-P + Pi + P2) + r) 

(2.113) 
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To conclude, in order to obtain approximations for the causal correlation functions 
C2, C3 and C^, we have first to solve the TDHB equations (|232|) and (p33|) for a^^\t) and 
and to evaluate the matrices r, R, t and T for the TDHB solutions. We then solve 
equations (ICTj) and (|]9g) for /(f ,t) and m{t" ,t) - l{t" ,t)a^^\t) {t" > t > to) backward 



in time from t" to to with the boundary conditions (|2.91|) - (|2.99|) . Finally, we evaluate 
expression ( |2.105|) . Therefore, it is not necessary to solve the dynamical equations (|2.79| ) 
and (ICTp for /?(t,f' ) and S(t,t") . 



The expressions we obtain in this way for three-point and four point correlation 
functions are not completly satisfactory. Indeed, in the static case, when a;^°^(t) = ao and 
H'^°^(t) = Ho are time-independent, the approximate correlation functions do not depend 
only on the time difference t' — t" as they should (the two-point function in the symmetric 
phase is a particular case). We will cure this artefact by optimizing with respect to the 
initial state in section 3. 

Retarded Green Functions : 

The function pf is an approximation for the Green function with the time-ordered 
Green function : 



C|^(x, y, t, t") = Tr (T$^(x, t', to) $^(y, t" , to)D{toj^ 

-Tr(<|.^(x,t',t")^(to)) Tr($^(y,t",to)I^(to^^ ' ^' ^ 



To obtain the approximations for the anticausal two-time functions, it is sufficient to 
replace in the functional ( |2.36| ) the term Kc by the analogous term Ki,. Formula ( p.l05| ) 
and (|2.106|) are then interchanged. 

The retarded Green function is defined by : 

X$<E.(x,y,t',t") = -tQ{t' -t")Tr {^,t', to), <^''{y,t",to)]D{to)) . (2.115) 

It is interesting to obtain an approximation for the retarded Green function. It is this 
Green function which appears in the study of the response of the system to a small exter- 
nal perturbation. Moreover, the KMS condition allows to relate the different components 
of the Green function in the real time formalism for Quantum Field Theories at finite 
temperature |r7|. These components can all be expressed from the retarded Green func- 
tion. The real parts of the time-ordered Green function and of the retarded Green 
function x coincide. They describe phenomena such as Debye screening in a plasma in the 
problem of charged particules in an external electromagnetic field. But one has to con- 
sider the imaginary part of the retarded Green function to describe damping phenomena 
|T8[| . The optimization with respect to the initial state in section 3 will allows to obtain 



approximated Green functions wich verify the KMS condition (the fiuctuation-disspation 
theorem). 
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From expression ( p.l06| ) and the analogous formula for the correlation function with 
the anti-T product, we obtain an approximation for the retarded Green function : 

xM^,y,t',t") = -^&{t'-t") [-4(m*(t",to)-«o/*r,to)) r (m*(t', to) - to)«o) 

+tr to)So/*(t', to)r) - tr to)r/*(t', to)Eo] ^ 



In the symetric phase : 

, y, t', t") = i e{t' - t") 4 m*(t", to)rm*(t', to) 



(2.116) 
(2.117) 



The exact expression (|2.115|) remains identical if we translate the refence time to 
to another time t such as t" > t > to. By using the dynamical equation for S(°)(t) 
and equations ( |2.89|) and (|2.90|) , we show that the approximation (|2.116|) has the same 
property. Indeed, if we call R{t',t") the expression obtained by replacing in the bracket 
of equation ( |2.116D to by t with to < t < t", we have for t' > t" : 



dR 
Itt 



2(m*(t',t) -aW(t))/*(t',t) 



5{t - 1") 



(2.118) 



We deduce : ^ = pour t <t" . We thus have 



X#4.(x, y, t', t") = -i 0(t' - t")R(t', t" 



(2.119) 



When the time t' tends to t", the exact formula ( p.ll5| ) coincides with the commuta- 
tor of two fields ^^(t', to) at equal time which is vanishing. The approximated expression 
gives also xi>*(t',t') = since mf(t",t") = and /*(t",t") = 0. 

Similarly, we can obtain approximations for the retarded Green function with two 
field operators defined according to : 

X**,** (x, y, z, u, t', t") = -id{t'-t")Tr ([<I>^(x, t', to)$^(y, t', to), $^(z, t", to)$^(u, t", to)]D(to)) . 

(2.120) 

In the symmetric phase, from expression ( p.ll3| ) and from the analogous formula for the 
correlation function with the anti-T product, we deduce an approximation for the retarded 
Green function with two field operators : 

X$$,<i.$(pi,P2,P3,-X1p»'^''^") = -^^(^' -^") / '^^'^ 

tr[/**(p - pi - P2, -p, Pi, P2, t", to) Ho(p) ^**(P, -p + Pi + P2, P3, - Pi, t', to) t] 

-tr[/**(p - Pi - P2, -P, Pl, P2, t", to) T /**(p, -p + Pi + P2, P3, - Pi' ^o) So(-p + Pi + P2)] 

(2.121) 

Again, the right member of this equation does not depend on the time to. Indded, if 
we call S{t', t") the expression obtained by replacing in the right member to by t, with 
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to < t <t'' we have for t' > t" : 



It 



and we deduce : ^ = for t <t" . We have therefore : 



So] 5{t - t") , 



X*#,4.$(Pl, P2, P3, - E P- = ^(^' - OS{t\ t") 



^2.122) 



(2.123) 



The limit when the two times t' and t" coincide : 



When t' ^ f', t' = f" = t + 0, expression ( |2.105| ) for the approximated two point 
function becomes : 



(3f{t,t) = 2(m*(t,to) - dor{t,to)) (Ho - r) (m*(t,to) - r(t,to)«o) 



+-tr 
2 



/*(t,to)(Ho-r)/*(t,to)(Ho + r) 



;2.124) 



At the initial time t = to, we get as we should the TDHB result : Pfito, to) = Go- Indeed, 
mf(to,^o) = and /*(to,^o) = 0. However for later times, this approximation is in 

general different from the result of the TDHB evolution : ^S^^^t) = G^^\t). The case at 
equilibrium in the symmetric phase is a special case since, as we will show in section 3, 
we have : 

/3f(p,At = 0) = Go(p) (2.125) 

This won't be true in the asymmetric phase. The result ( p.l25| ) is to relate to the fact 
that, at equilibrium in the symmetric phase, the two approximations for the two-time and 
two-point function obtained with or without the optimization with respect to the initial 
state coincide (see below eq. ( |3.58| )). 

Let us examine the four-point function when the two times coincide. In the sym- 



metric phase, the limit t' = t" — = t + in expression ( |2.111| ) for the four-point function 
Sff gives : 

Sff (pi,P2,P3,-EPi,^,^) = J ^^^^^[^**(P-Pi -P2,-p,Pi,P2,t,io)(Ho(p) -r) 

^**(P, -P + Pi + P2, P3, - E P»' ^0, t)(Ho(-p + pi + P2) + t) 

(2.126) 

The boundary conditions for the function write 



^n*(Pi,P2,P3,-EPi'^'^) 



<5'(Pl + P3) +5'(P2 + P3) 



(2.127) 



the other matrix elements of Z** being equal to zero. We thus have at the initial time : 



Sn*(pi,P2,P3,-EP«'^o,to) = - H?i(-pi)H?i(p2) 



<^^(Pl + P3) + 5^(P2 + P3) 



(2.128) 
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The mean-field approximation for which C4 = corresponds to (see eq. ( |2.74| )) : 

Sff (x, y, z, u, t, t) = 2/3f (x, z, t, t)/3f (y, u, t, t) + 2/3f (x, u, t, t)/5f (y, z, t, t) 



(2.129) 



= 2GW (x, z, (y , u, t) + 2GW (x, u, t)^^ (y , z, t) 
In momentum space, this corresponds to : 

Sff (Pl, P2, P3, - E = 2/5f (-P1, t, t)/?f (P2, t) [5(P3 + Pi) + 5(P3 + P2)] (2.130) 



At the initial time, we recover the result of the Wick theorem given by expression p.l28| in 
momentum space. However, for later times, owing to the correlations we have taken into 



account in the dynamics, approximation ( |2.126|) for the four-point correlation function is 
different from the result of the Wick theorem. 

In nonrelativistic many-body problems, the mean-field TDHF calculations correctly 
predict the expectation values but underestimate the fiuctuations. Balian and Veneroni 
have showed the advantage to use a variational principle adapted to calculate the fiuctua- 
tions. For a pur state in the case of heavy-ions collisions, one obtains in this way a better 



agreement with experiments for the width of the mass distributions |19]. The example of 
the quartic oscillator in one dimension shows also the usefullness of the method developped 
by Balian and Veneroni to calculate the fiuctuations of the operator x^, < > — < >^, 
though staying in a mean-field approach [|14|. 
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3 Optimization of the initial state 



3.1 Introduction 

In physical situations, the initial state is not completly known. Only the expectation 
values of a few observables are known, for instance : 

(<l>(x)),,^ , ($(x)<|.(y)),^ , ($(x)$(y)$(z)),^ , ($(x)$(y)$(z)$(u)),^ . (3.1) 

The exact state D{to) can be written as exp{—i3H). The operator H is equal to 
the Hamiltonian of the system at thermal equilibrium but it can also contain other terms 
corresponding to conserved quantities, like the charge operator for a complex scalar field. 
For a nonequilibrium situation, the operator H is different from the Hamiltonian H. 

We will approximate this initial state in the following way. We use the fact that the 
operator e~^^ is a solution of the Bloch equation : 

A e""^ = -H e""^ , 0<u<p . (3.2) 
du 



Following the general method of 0, [|T^ to evaluate variationally a given quantity, we 
introduce an ansatz operator T>{t) and we consider equation (|3.2|) as a constraint for the 
complex times t = tQ + i{i3 — u) with the boundary condition V(tQ + i(3) = 1. We introduce 
also a Lagrange multiplier A{t) associated to this constraint. The functional to minimize 
writes : 

Z{Ait),V{t)) = Tr{Aito + tO)V{to + tO))-Tr dtAit) + t H V{t)\+Zdyn , 

(3.3) 

where Z^yn is given by expression ( |2.15|) of section 2. The ansatz operators P(t) and A(t) 



are subject to the boundary conditions ^(oo) = 1 and P(to + i (3) = 1. 

Let us note that for a complex value of t, A{t) is not close to 1 in the limit of vanish- 
ing sources. Indeed, the stationarity condition of Z with respect to arbitrary variations 
of T^lt) writes : 

- iA{t)H = (3.4) 

and the solution of this equation for t = + i{P — u) is equal to : A(t) = A(to) exp(— (/? — 
u)H). 

The approximate dynamical equation for T>it) and A(t) are coupled. V^to) depends 
on A{t) for t > to and therefore on the sources J contrary to the exact state D{to). 
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3.2 The variational spaces 



For the variational spaces for T>(t) and A(t), we choose again Gaussian operators charac- 
terized by the contractions ad(x, t), Hd(x, y, t) and aa(x, t), Sa(x, y, t) where the time t is 
on the two segments [to + i P,to[ and [to, +oo[ in the complex plane. 

The boundary condition P(to + iP) = 1 writes : 

1 



(3.5) 



Let us write the functional Z (|3.3|) with the elementary contractions : 



rto 



f 

Z{na, aa, S^; Ud, ad, H^) = n(to) - dt 



or, by integrating by part : 



dn{t) 
dt 



A{t)=cte + i n{t) / d^x £c{^, t) 



(3.6) 



f 

Z{na, aa, Sa; Ud, ad, H^) = n(to+iP)+ / dt 



dn{t) 
~ltt 



(t)=cte - i n{t) I d^x Scix, t) 



We have used the notation : 

J d^x ^,(x, t) = Tr {V{t) A{t)H) =< H >, 



~\~Zdyn 

(3.7) 
(3.8) 



3.3 The dynamical equations 

For complex values of the time to + i P < t < to, the dynamical equations for the contrac- 
tions ad, '^d, aa and write : 



2idd = - (Hrf + r) He (H^ - r) r a^.c + (H^ + r) 

iEd = -{Ed + r) Tic (2d - r) , 
2i aa = (Ha - r) Tic (2a + r) r a^-b - (Ha - r) , 

i Sa = (Ha - r) He (Ha + t) . 



(3.9) 

(3.10) 
(3.11) 
(3.12) 



These dynamical equations are coupled through the boundary conditions V{to + il3) 
l,A{t = +oo) = 1 which write : 



H^^(to + ^P)=0 , —adito + if3) = Q 



E-\t = +oo) = , ^aa(t = +oo) = . 



(3.13) 
(3.14) 
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For the contractions a^, et etc, ^c, the dynamical equation write : 



db = , (3.15) 

Hfe = , (3.16) 

i ac = T Wc , (3-17) 

tE, = 2 (s,Kr-rKHc) . (3.18) 



We have '^^^^'^ = 0. We have introduce the vector Wc and the matrix Tic (analogous 
to w and Ti ) to express the variation 6 < H > in terms of the variations 6ac and (5Hc. All 
the contractions a et S are continous at t = to- In the complex time interval [to + ip, to], 
ab and Eb are constant : = ab(to), = Eb{to). 

3.4 Expansion in powers of the sources 

We use the index ^'^^ for the solutions of the dynamical equations in the limit with no source 
: a,{t) = af\t),E,{t) = Hj)(t),a,(t) = a^^\t),Ea{t) = S(o)(t). Since Ait) ^ 1 even in 
the limit with no source, this limit does not correspond to 4-<^a = 0, E'^^ = 0, contrary 
to what happens when we do not optimize the initial state (see section 2). Therefore, 
Ef^ 7^ E^^^ and a^^ ^ af^ . However, we can use a continuity property at t = to- in the 
limit with no source, A{to + ^0) = ^(to + 0) ^ 1. At the time to, we thus have : 

af{to) = af{to) = af{to) = ao, (3.19) 

Et\to) = E^\to) = Ef{to) = Eo, (3.20) 
where oio and Sq are the HB static solutions associated to the operator II and verify : 

w;o = , (3.21) 

So Ho r = r 7^0 So . (3.22) 

So et 7^0 are related by equations similar to ( p.55| ) and ( |2.56|) . The expressions for 
the vector wq = w(q;o,So) and the matrix Hq = 'H{ao,Eo) are obtained from those in 
Appendix D by replacing /o and go by the time-independent functions fo and ^o- (We 
recall that Ti. comes from the variation oi < H > with respect to S and has to be evaluated 
at tto. So to obtain Ho- Since H is a time-independent operator. Ho is time-independent.) 
We stress the difference between the quantities S^°\ Sq and So. 
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The expansion at zero order of equations (|3.9|)-( |3TT^ ) writes : 

2^aT = {^T+r) , (3.23) 



2 ^ = - r) , (3.25) 
zS(°) = -(si°)-r) 7^o (s(°) + r) . (3.26) 



Since = 0, we have for to + il3 < t < : af^ = and = 0. By using the 
continuity at t = to, we thus obtain for < u < (3 : 

a^d\'to + i(3 — i u) = . (3.27) 

{^aa)''^\tQ + i(3-iu) = Q , (3.28) 



The solution of ( p.24|) with the boundary condition {tQ+ijS) = is the following 



2^°^ (to + 1(5 -iu) = -T coth(M Ho r) , (3.29) 
which can also be written as : 

S?)(to + z/3 - zn) = r f 1 - f ) (3.30) 

Y 1 — exp(— n7to2r) J 

or in the uniform case : 

row s /-coth(^) \ , , 

\ -co'pcoth(^) / 



with, for the theorie, ujp = y—goip) and 

^o(p) = -[p' + ml+ ^-ifl + ^Go(x, x) j (3.32) 

At u = 13, we have to obtain the static HB solution associated to the Hamiltonian H , Sq 

Ef^ {to + i 0) = Ho = -r coth(/5 Hq t) , (3.33) 



The solution of ( |3.26| ) writes : 



S^°)(to + ip-iu) = -T coth((/3 - u) Ho r) . (3.34) 
One checks that for u = (3, E!-^^ ^(to) = 0. 
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At the first order, w^^\to) = wo(to) = and E^^\to) = Hq. Therefore the right 
members of equations ( p.l7|) and (|3.18|) are vanishing at to- We deduce that at the first 
order in the expansion in powers of the sources, ac and Sc are constant in the interval [to + 
iP, to] and equal to ao and Sq respectively. We can also obtain this result from the relations 
of Appendix C and from the expressions obtained above for aW(t), Ef\t), E^^\t). 



At the next order, the approximated initial state characterized by 0:^(^0), ^^(to) is 
no more equal to the static HB solution ao, So. 

Instead of using the expansion of equations (|3.9|) and (|3.1CI|) for ad and S^, it is 
more convenient to expand at the first order equations (|3.17|) and ( p.l8|) for ac and Sc. 

For to + i(3 < t < to : 

1 



^^^f3'^{x,y,tX: 



2 ) 



dt ^ 



2 

S* Tio r - r Tio S* 



-— ^0 R 
2 



[kl ^fk) T 



IrR 

2 



{kl ^/fe)-0 



(3.35) 



(3.36) 



The matrices t, T and f, R associated to the operator H and depending on So and on 
ao are analogous to the matrices t,T,r and R. The quantities Ho, ao,So are time- 
independent for to + iP <t <to- The matrices i, T and f, R are also time-independent for 
to + iP < t < to- Indeed, these are defined from the expansion of Wc and He in terms of the 
variations 6ac and 5Sc (see equations (|2.77| ) and ( p.78| )) and they have to be evaluated for 
the values a^^-* and S^^-*. Now, as we have noticed before, these quantities do not depend 
on the complex time and are equal respectively to ao and Sq. For the same reason, Sq 
appears in the last two lines of ( |3.36|) and not S^°^ as in equation ( p. 801) which describes 
the evolution of S*(t,t") for real time. 



3.5 Approximation for the two-point function 



In the following we consider an uniform problem. In momentum space, equation ( 3.35| ) 
for writes : 

^ ft*(P. t, n = 5o(-P) ,<?f (P. t, n + 1 ^. ll sf.(P.. P.. P. t, n . 

(3.37) 

t^^Pfip,t,t") = -p!{p,t,t") (3.38) 
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Symmetric phase 

In the symmetric phase, i^q = and the equation for the function /3f(p, for 
to + iP < t < to writes : 

^Pf(p,t,t") = goip)PfiP,t,t") , (3.39) 

where 

- ^o(p) = p^ + ml + -Go(x, x) (3.40) 
We can therefore relate /3i*2(P> ^o, t") to /3*2(P> + ip, t") : 

(^l) (p,to,n = ^(P) (j^l) (p,to + ^An (3.41) 
where ^(p) is the matrix : 

, / coshu;p/? \ , , 

^(P) = . / . J (3.42) 
ycUpSmhcUpP coshcjpp / 

with Up = \/-(7o(p). 

We then use the fact that T'(to + ^/3) = 1. We deduce : 

adto + iP) = Mto + if3) , ^c{to + iP) = ^b{to + iP) ■ (3.43) 

By introducing the functions 6 and A which come from the expansion at the first order in 
powers of the sources of at and Sf, (and which are analogous to /5 and S), we thus have : 

p{to + 1(3, t") = 6{to + 1(3, t") (3.44) 

S(to + i(3, t") = A{to + i(3, t") (3.45) 

Now ab and are constants for to + i(3 < t < to (eq. ( |3.15| ) and ( p.l6| )). Consequently : 

5{to,t") = 5{to + i(3,t") (3.46) 

A(to,t") = A(to + ^At") (3.47) 

Let us compare the expansion of H^, He and at to- By using equation ( |C.6| ) at t = to, 

Afc(to,t") = -(So + r) kit", to) (Ho - r) , (3.48) 
and from equation ( |C.lj ) : 

6kito,t") = {% + T){mk{t",to)-lk{t",to)ao) ■ (3.49) 
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The comparison of these last two equations to ( |2.86| ) and ( p. 87] ) leads to : 

h{to, t") = Pk{to, t") + 2r(mfc(t", to) - to)«o) (3.50) 

Afc(to, t") = J:k{to, t") + 2 %lk{t\ to)T - 2 rhit", to)So . (3.51) 

Finally, we use : 

(3{to + if3,t") = 6{to,t") (3.52) 

j:{to + if3,t") = A{to,t") (3.53) 
In the symmetric phase, /* = 0, therefore : 

/?*(to + t") = /3*(to, t") + 2rm'^(t", to) (3.54) 

From equation ( 3.41|) , we obtain : 

/3*(p,to,t") = {A-\p) - l)"'2rm*(p,t",to) (3.55) 



For t" > t' , the relation ( |2.109|) of section 2 is therefore modified according to : 



/?f (p, t', t") = 4 m*(p, t', to){A-\p) - 1)- Vm*(p, t" , to) (3.56) 

with 



2 coth(|a;p) ^ 



(•^-(P) - 1)-' ^ ( -^T^^'-^M (3.57) 



For t' > t", we have : 

/?f (p, t', t") = 4 m*(p, t', to) ((A-i(p) - l)-i + l) rm*(p, t", to) (3.58) 

We check again that in the limit t' = t" = to, we obtain the TDHB solution at to : 
/3f(p,to,to) = ^ coth{f3ujp) = Go(p,to). 

For the two-point function in the symmetric phase, one can check from expression 
( |2.57D for the matrix Sq (and where the quantities which appear correspond now to H 
and no more H) that : 

Ho(j9)-r = 2(A-i(p)-l)-V (3.59) 
Therefore expressions (|2.109|) and (p.56|) for f3f in the symmetric phase obtained without 



or with optimization with respect to the initial state coincide when Sq = Ho which is 
the case at equilibrium. This is not true in the asymmetric phase. It will be also more 
interesting to study the four-point function in section 3.6. 

It is also interesting to calculate the approximation obtained for the retarded Green 
function ( |2.115 ) when we optimize both the dynamics and the initial state. In the sym- 



metric phase, we have for the time-ordered correlation function for t' > t" : 

/3f(t',t") = 2m*(t",to)/5*(to,t') (3.60) 



33 



For the anti-T product and t' > t" : 

/3f(t',t") = 2m*(t',to)/5*(to,t") (3.61) 
By using ( p.55|) and the fact that 



{A-\p) - + [t{A-\p) - = -T , (3.62) 

we obtain : 

X**(p, t', t") = i6{t' — t") 4m* (p, t", to)Tm'^{p, t', to) (3.63) 

which is identical to the formula ( p.llTp obtained without the optimization of the initial 
state. For the retarded Green function the approach of section 2 is therefore sufficient. 

The case at equilibrium 

In this case, the two-time functions should depend only on the difference t' — t" . 
Unfortunatly, this is not true for the expressions we gave for the three-point and four- 
point functions in section 2. At equilibrium, a^^\t) = aQ,'^^^\t) = Sq and the matrices 
t, T, r, R are constant. In equations (p.89|) and (|2.90|) for and m^, the shift t' ^ t' + St is 
therefore equivalent to a shift of the initial time to to — 6t. By writing the corresponding 
shifts of the functions I'' and m'^ — I'^a^ and by inserting them in eq. ( ^.105| ) or in eq. 
( p.lOTp , one sees that, in general, our approximations are not invariant under the shift 



t' ^ t' + 6t, t" — > t" + 6t, to being constant. They contain a spurious dependence on the 
initial time to. 

The case of the two-point function in the symmetric phase is special. From equation 
( p.90|) for m*, we have (5m* = itrm'^dt and the variation of expression (|2.109| ) is equal to 



zero because 5(3f is proportional to —Sjj^ (70 + 222 which is vanishing for the static solution. 
However, this is no more true for the four-point function E** given by equation ( |2.113|) . 
The variation (5S involves two types of terms : one with the matrix R proportional to 
the coupling constant and another with the matrix TYq. Generally, 5S is non vanishing. 
This spurious dependence on the initial time will disappear when we take into account the 
optimization with respect to the initial state : the two-time correlation functions depend 
only on the time difference t' — t" as they should. 

At equilibrium, in the symmetric phase, we can solve equation ( p.90| ) for m* : 

mf (p, t', to) = cos(u;p(to - t')) (3.64) 

mf (p, t', to) = ^ sm{ujp{to - t')) (3.65) 



We define At = t' - t" . For At < 0, equation ( |3.56| ) gives 



/?f (p. At) = —(coth ^ cosu;pAt - - sinu;pAt) (3.66) 



2ujp 2 ^ i 
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For At > 0, we use equation ( p.58| ) : 

/?f(p,At) = -— (coth^ coscUpAtH-- sincUpAt) (3.67) 

We check that /5f depends only on /St = t' — t" . Expressions (|3.66| ) and ( |3.67|) are 
interchanged by changing At to —At, which was not exphcit on formula ( |3.56D and ( |3.58| ). 

For At = 0, we obtain /3f (p, At = 0) = G'o(p)- In this case, at equal time, the two- 
time and two-point function coincides with the mean-field HB solution for the two-point 
function. 

By taking the Fourrier transform in time of equations ( |3.66| ) and ( p.67| ) : 

/O r+oo 
d(At)e'(--^^)^*/3f (p. At) + / d(At)e^(^+^^)^*/3f (p. At) , (3.68) 
-oo JO 



we obtain : 



/3f(p,^) 



rip -Fl rip 



uj"^ — {ujp — iey uj"^ — (cjp -|- ieY 



(3.69) 



where n^, is the occupation number : 2rap + 1 = coth 

We will now show that our approximation for the two-time correlation function with 
two field operators satisfies to the fluctuation-dissipation theorem. This theorem can be 
written as a relation between the Fourrier transform in time of the correlation function 
C|^(x, x', t, t') and the response function 

X2(x,x',t,t') = rr([$^(x,t,to),$^(x',t',to)]/^(to)) (3.70) 

in the following way : 

X2(x,x',cu) = (l-e-^'^)C'|^(x,x',u;) , (3.71) 
where the Fourrier transform in time is defined according to : 

rf(t_t')eMt-t')^^(x,x',t-t') (3.72) 

-oo 

The correlation function gives us informations about the fluctuations in the equilib- 
rium state while the response function contains informations on the dynamics of the 
system which has been driven from the equilibrium by a small external perturbation. 
The function appears also in the dissipation energy of the system. Therefore, the 
fiuctuation-disspation theorem tells us that, in systems near equilibrium, the transport 
properties which are linear in the external forces can be calculated from the fluctuations 
at equilibrium [pOf . 

By using the definition 

Tr($^(x, t', to)$^(y, t", to)I^(to)) - y.(x, t')^{y, t") <|.(x, r)<|.(y, r) > (3.73) 
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with T = t' — t" , the fluctuation-dissipation theorem can be written in the following way : 

/+00 
dre"^^ < <l>(y)<l>(x,r) > 



hoo 



(3.74) 

dr e^'^(^-*^) < $(x, r - i/?)$(y) > 

In order to check the fluctuation-dissipation theorem on our approximations for the two- 
point correlation function, it is therefore sufficient to check that expression ( |3.66|) for 
At = r is equivalent to expression ( |3.67| ) for At = — r — ijS, which can be done easely. 



At equilibrium, expression (|2.117|) for the retarded Green function with two field 



operators in the symmetric phase gives : 



X<s>^{p,t',t") = -9{t' -t")— sin ujp{t' -t") (3.75) 



UJp 



3.6 Four-point function : a perturbative expansion 

Still in the symmetric phase and for the uniform case, we consider now the four-point 
function S**. 



3.6.1 The response function Ilji{q^ ,t' ,t") 

When one couples a source J**(x, y,t') to the the operator $^(x, t', to)'^'^(y, i", ^o); the 
linear response theory gives the variation of (^(x, y, t') = Tr($^(x, t', tQ)^^{y, t', to)D(to)) 
in presence of the source : 

/'4-00 r 

6G{x, y, t') = / dt" / dh d^u *-i.(x, y, z, u, t', t") J**(z, u, t") (3.76) 

Jto J 

where x**,** is the retarded Green function with four field operators defined in eq. 
( p.l2C ). In presence of a source of the following type : 



J**(xi,X2,t") = J^''{t")e-'^ •^^5=^(xi - X2) , (3.77) 

n+OO 
'to 

By introducing the Fourrier transform , we obtain for x = y : 



/'+00 r 

SG{x, y, t') = / rft" J**(t") / d^z e-'"^ -"x** *-i.(x, y, z, z, t\ t") (3.78) 

Jin J 



5G(x,x,t') = e-^^-/ dt"J**(t") / ^^;^^^,^^(k,-q-k,pi,q-pi,t',t") (3.79) 
The polarization function Ili^ is defined according to : 

(3.80) 



2 ,/ ,//^ ^ 6G{^,^,t') 



J**=0 
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1. e. : 



(27r)3 (27rj3 ~^ ^ P^' ^ ~ P^' ^" 



(3.81) 



A variational approximation for the polarization function Tin is obtained by using expres- 
sion ( |2.121[ ) for x-ixJ-,** ^ function of and Sq. 

In the static case, n^(q^, t', t") depends only on At = t' — t" and is proportional to 
9{t' — t"). We thus define its Fourrier transform in time by : 

^^(q^ u;) = / d{At) e^-^*+^'^*^H(q^ At) (3.82) 
Jo 

The polarization function n/j(q^, t', t") can be recognized as the Fourrier transform 
with respect to x — y of the retarded Green function for the operators $^(x) and $^(y) : 

X$^^2(x, X, y, y, t\ t") = -ie{t'-t")Tr ([$^(x, t\ to)$''(x, t', to), $''(y, t", to)$''(y, t", to)]/^(to; 

(3.83) 

We have : 



(3.84) 



^3^e-q.(x-y) ;^^,^^,(x,x,y,y,t',t") 

= n«(q',t',n 

This function is directly related to the viscuosity in the "l*^ theory. From the spectral 
representation : 

X$2$2(p,po) = lim / duj — , (3.85) 

e^0+ J-oo Pq — UJ + le 

the viscuosity is defined by : 



r]q,2 $2 = lim 



P<i.2,4.2(p,po) 



(3.86) 



Calculations of rii^2 q,2 have been done in the limit of high temperature and the result have 



been compared to the classical approximation pT 



Actually we will begin by calculating the following non-retarded Green function : 

(3.87) 



n(q2,t',t")^L^^ii(''''''^') 



J**=0 



which can be written with the four-point function Sf* as : 



^(q^t^t"^ 



S'k f d^pi 



Sf*fk, -q-k,pi,q-pi,t',t" 



(3.^ 



(27r)3y (27r)3 " 

In the following, we will therefore consider the function S**(k, — q — k, pi, q — pi, t', t"). 
We will show the differences between the approximations for the function n(q^,t',t") 
when we do or do not optimize with respect to the initial state. 
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3.6.2 Resolution of the dynamical equations for complex values of the time 



In the symmetric phase, we have for + i(3 < t < to : 

d 



z^Sff (pi, p2, P3, - E P- n = -^t2 - ^li (3.89) 

^|sf* = ^7o(-P2)Sff - - / (3.90) 

= Upi)^u - - l^'iM) I (3.91) 

^|sf2* = 9o{-P2)^ti + 9o{Pi)^t2 (3.92) 
where I is the following integral : 

i=[ ^?i*(p4.P5,P3,-Ep*'*'^") (3-93) 

"'P4+P5=pi+P2 {27T)'^ [Zny ^ 

These equations are analogous to ( |2.80| ), the matrices Hq, r and R being replaced by the 
time-independent matrices Hq, f and R associated to H and E^^^ by Sg. We see that we 
are not abble to obtain directly an expression analogous to ( ^■41[ ) which would relate the 
matrices S(to,t") and S(to +iP,t") since the loop term |S5i / involve all the momenta. 
We will use a perturbative expansion where we neglect this loop term to lowest order. 
From the expression for S^i, we see that this expansion will be valid when 

^coth^«l (3.94) 
where m(/3) is the self-consistent mass defined by : 



At zero temperature, the condition (|3.94|) for the validity of the perturbative expansion 



writes simply | <^ 1. At high temperature it becomes : | <^ /5m(/5) with m^(/?) ~ 

III' ^ 24/32 . 

To lowest order, by neglecting the loop term I, the solutions of equations ( 3.89| )- 
:92D write : 



Ef *(k, — q— k, pi, q— pi, M, t") = Ai cosh(i73-u)+i?i sinh(f23M)+Ci cos\i{Vl/^u)+Di sinh(n4M) 

(3.96) 

Ef2*(M, t") = cjq+k [^1 sinh(n3M) + Bi cosh(fi3M) + Ci sinh(il4-u) + Di cosh(i74'u)] (3.97) 
S|*(u, t") = tUk [Ai sinh(r23M) + Bi cosh(r23M) — Ci sinh(r24n) — Di cosh(r24M)] (3.98) 
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1^22 (u, t") = cUkCUq+k [^1 cosh(n3M) + Bi smh{il3u) — C\ cosh(f24M) — D\ sinh(f24'u)] 

(3.99) 

where we have used the variable u defined by : t = tQ\%{fi — u) and we have introduced 
the frequencies 

^3=^^q+k+^k (3.100) 
fi4 = '^q+k-'^k (3.101) 



cuk being defined by cuk = y ^5'o(k). 

We can therefore express E**(k, — q — k, pi,q — pi,to + 'iP^t") as a function of 
S**(k, — q — k, pi, q — pi, to, t"). The result is the following : 



S**(k, -q-k, pi, q-pi, t") = A-i(k)S**(k, -q-k, pi, q-pi, to, t") {A' )-\q+k) 

(3.102) 

where the matrix A{k) is defined by ( |3.42|) : 

A{k)- 



cosh(^k/?) ^i^^ 
ciJk sinh(ti;k/5) cosh(co'k/5) 



A-\k) 



cosh(c<Jk/3) 
-cuk sinh(ti;k/5) cosh(co'k/5) 



sinh{ciJi5./3) 



(3.103) 



(3.104) 



By introducing a column formed with the elements Sjj, we can write this expression 
in the following way, which will be usefuU in section 3.6.4 : 



(k, -q-k, Pi, q-pi, to, t") = A^(k, q+k) 



/ ^11 \ 

^12 
^21 



(k, -q-k, pi, q-pi, to+i/3, t") 

(3.105) 



where A^(k, q + k) is the following 4 by 4 matrice : 



c 



B 



iV(k,q + k) 



_D 

D 



\ 



^k 

c 



with 



' <^k+q 
^k+qC A 

UkD -^B A 

VcUkCUk+qfi LOl^D C^k+qC A / 



A = cosh(co'k+q/3) cosh(u;k/5) 
B = sinh(co'k+q/?) sinh(ci;k/3) 
C = sinh(c<Jk+q/3) cosh(co'k/?) 
D = cosh(ci;k+q/?) sinh(tt;k/5) 



(3.106) 



(3.107) 
(3.108) 
(3.109) 
(3.110) 
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If we do not optimize with respect to the initial state, an approximation for the 
four-point function T,f^(t',t") is given by equation ( p.lllj ) for t" > t' . 

On the other hand, if we optimize with respect to the initial state, Ef*(t',t") is 
given by ( gTTO| ) for t" > t' : 



S?f(Pi,P2,P3,-X!P»'^''^") = - / 



tr 



X 



/**(p - pi - P2, -p,Pi,P2,t',to)S**(p, -p + pi + P2,P3, -XlPi'^O,^") 



(3.111) 



and we use ( |3.51|) and ( p.53| ) to write 

S**(P1, P2, P3, - P»' ^0 + t") = S**(pi, P2, P3, - Pi' ^0, t") 



+2 



2o(-Pl)^**(Pl, P2, P3, - Y Pi' ^"' ^o)^ - ^^**(Pl, P2, P3, - Y Pi' ^"' ^o)2o(P2) 

(3.1 



2) 



or : 



A-i(k)E**(k, -q - k, pi, q - pi, to, t")(^^)-'(q + k) - S**(k, -q - k, pi, q - pi, to, t") 
2 [So(-k)/**(k, -q - k, pi, q - pi,t", to)T - r/**(k, -q - k, pi, q - pi, t", to)So(-q - k) 

(3.113) 

where So is given by ( p. 57] ), go being replaced by ^o- 

The matrix elements S**(to,t") can therefore be written in terms of /**(t",to)- By 
using ( p.lll| ), we finally obtain an approximation for the two-time and four-point function 
Ef*(t',t") which involves only the function We recall that this approximation has 
been obtained by neglecting the loop term I ( ^.93| ) in the evolution equations of S** for 
complex values of the time. 



3.6.3 Resolution of the dynamical equations for / in the static case 

The dynamical equations for the matrix elements of / are given by (|]89D. In the static 
case, E^^\t) is equal to So given by ( p.57|) , a^^\t) = ao and go is time-independent. In 



the symmetric phase , equations (|2.89|) write more explicitly in momentum space : 

^^/f *(k, -q - k, p„ -p, + q, t", t) = -<7o(k)/2? - ^?o(q + k)/** + J (3.114) 

^jfi2=lti-9Mlt2 (3.115) 
,|/|* = /f*_^„(q + k)4<^ (3.116) 

^J^lt! = lt2+ltl (3.117) 
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with the boundary conditions ( |2.127| ). J is the following integral : 

= -! / (0 [So(p)/**(-p, -q + P, P., -P. + q, t", t)r 
-rZ**(-q - p, p, Pi, -pi + q, t'\ t)Eo{p) 



(3.118) 



J 



A J (27r)3cjj 



coth(^) /f2*(-p, -q + p, Pi, -Pi + q, t", t) + ^2i*(-q - p, p, p*, -p* + q> t", t) 

(3.119) 

In the previous equations, in order to make the notations more compact, we have ommitted 
when it was possible the arguments of Let us remark that, on the contrary to 

equations (|3.89|) - (|3.92|) for S**, the solution Ho(p) appears inside the loop integral. 

Similarly to the resolution of the equations for S**, we will solve the equations 
for perturbatively, neglecting first the loop term J. In the following section, we will 
take into account of the loop term J to first order. To lowest order, the solutions of 
(|]Tl|)-(|]TT3) write : 

^n*(k, -q-k,pi,q-pi,t',t) = - 



li2 (k, -q-k, pi , q-pi ,t',t) 
^21* (k, -q-k, pi, q-pi, it', t) 
1^2 (k, -q-k, pi , q-pi ,t',t) 



1 



1 



6%k + pi) + 5^(-q - k + pi)] cos(cJk+q(t-t')) cos(cJk(t-t')) 

(3.120) 

5^{k + pi) + 5\-q - k + pi)] sin(cjk+q(t-t')) cos(cJk(t-t')) 

(3.121) 

5^(k + pi) + 5^(-q - k + pi)] cos(cUk+q(t-t')) sin(cuk(t-t')) 

(3.122) 

6\k + pi) + - k + pi)l sin(cuk+q(t-t')) sin(cuk(t-t')) 



(3.123) 

Formula ( ^3Tl| ) and (|3lT3|) then give, for At = t' - t" < : 

Eff (k, -q-k, pi, q-pi, t', t") = 2 [^^(k + p^) + <53(-q - k + pi)] /5f (k, At)/?f (q+k. At) 

(3.124) 

where jSf (k. At) is given by ( |3.66| ) . We thus have at the lowest order : 



nW(q2,At) = 2z J ^/5f(k,At)/3f(q + k,At) 



(3.125) 



When we neglect the loop terms I and J in the dynamical equations for S** for the 
complex time and in the dynamical equations for for the time t > to, we thus recover 
the Wick theorem as we should. Let us stress that we have neglected the interactions 
which appear in the terms I and J but the interactions remain in the frequencies Uk which 
involve the renormalized self-consistent mass m{f3) (eq. ( |2.66| )). 

At this order, Sf* depends only on the time difference At = t' — t" and our ap- 
proximation satisfies the fluctuation-dissipation theorem for the two-time and four-point 
correlation function. 
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If we apply formula obtained without optimization with respect to the initial 

state, with the solutions given by ( p.l2C| )-( 3.123| ), we find again the result of the Wick 
theorem. When we neglect the terms with the interaction in the dynamical equations for 
S** and the two approximations for Ef t") obtained with or without optimization 
with respect to the initial state therefore coincide. Accordingly, at zero order, we can use 
instead of expression ( p. 1131) the simpler expression : 



S (k, -q-k,pi, -pi+q,to,t") 



/**(k, -q-k. Pi, q-pi, t", to) (2o(-q - k) 

(3.126) 

The existence of such an expression with matrices (to distinguish from a vectorial expres- 
sion like (|3.105|) was not obvious at the beginning. 

For the polarization function n^(q^, t', t"), we obtain to lowest order : 



^g)(q^^^^")= -e{t'-t") J 



1 



(27r)3 t^pCUq+p 

(3uJr) A • A f^^O+q A • A 

coth cos Up At sm Wp+q At + coth — ^— ^ cos cjp+q At sm At 



(3.127) 



Its Fourrier transform in time writes 



[n 



p+q 



+ 1) 



q+p 



p+q 



{uj + iriY - (cjp+q + ujp) 



p+q 



{U + i7]Y - (C^p+q - UpY 

(3.128) 

A similar expression appears for instance in |^. Here, we have derived it variationally. 
We stress again that this formula is nonperturbative in the sense that it contains the 
renormalized self-consistent mass. 



3.6.4 Resolution of the dynamical equations with the interaction terms I and 
J at first order 

At first order, we will solve the dynamical equations (|3.89| )- p.92| ) for 
E**(k, — q — k, pi,— pi + q, M,t") by replacing in the loop integral I the solution E** 
obtained at the lowest order. We thus have to solve a first order differential sytem with 
a nonhomogeneous term which involves : 

with (see the analogous definitions ( p.lOO|) and (|3.101|) ) : 



r cP'l 

I{u) = / - — -r [a coshifl^u) + b smh.{Q^u) + ccosh.{QQu) + d smh.(flQu)] (3.129) 
J (27r)'^ 



^5 = ^q+1 + (3.130) 
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fig = — Oi!\ (3.131) 

The constants a,b,c and d depend on the momenta 1 and q. They can be expressed in terms 
of the matrix elements — q — 1, Pi, q— Pi, ^o, t") and cosh^fl^P), sinh^flQ j3) , cosh.{flQP), 

smhlfl^P). At the order we consider, we can use ( p.l26|) and express a,b,c,d with the 
matrix elements — q — 1, Pi,q — Pi,t",to), with coth(^^), coth(^^) and with the 

occupation numbers ni et ni+q defined by : 



2ni + 1 = coth(^) . 



(3.132) 



The relation between E**(to,^") and E**(to + iP,t") writes now : 



v-i$<J> 



(k, -q-k, pi, q-pi, to, t") = N{k, q+k) 



/ ^11 \ 

^12 

v-i$<J> 

^21 



(k, -q-k, pi, q-pi, to+^/?, t")+W{to., 



(3.133) 

where the matrix A^(k, q + k) given by ( p.l06| ) does not contain explicitly the coupling 
constant and W is a vector whose components are proportional to the coupling constant 
b. 



By using the fact that the first term can be written as ( p.l26| ), we write in matrix 
notation : 



S**(k,-q-k,pi,q-pi,to,n 



E:o(k)-r / (k, -q-k,pi,q-pi,t",to 



(So(-q-k) + r) +X{to,t") 



where 



and 



X{to,t") 



Xi X2 
X3 X4 



(3.134) 
(3.135) 

(3.136) 



X, = (1 - N)r/W, 

X(to, t") depends on the momenta k, k + q. The elements Xi can be expressed relatively 
simply with the integrals In, Jn, Kn, Ln, Mn of Appendix E : 



{ (^k+q + + 1) [-^^3(/i + I2 + Ki + K2 + Ks + Ki) 

- coth(^)(Ji + J2 - Ml - M2 - M3 - M4) 
(rik - rik+q) [-^^4(^3 + /4 + -^vi + L2 + L3 + L4) 

- coth(^)(J4 + J^-N^-N2-N^- 



} (3.137) 
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^2 = 8;^{ K+q + nk + 1) 



-^3 coth(^)(/i + /2 - - - 7^3 - K,) 



-{Ji + J2 + Mi + M2 + M3 + M4)] 
-(J4 + J3 + iVi + + A/'s + A^4 



fi4COth(^)(/3 + /4 - Li - L2 - L3 - L4) 



-{ ('^k+q+^k+l) 



-^3 coth(^)(/i + I,-K^-K2-Ks~ K^) 



(3.138) 



-(^1 + J2 + Ml + M2 + M3 + M4)] 
-(rik-r^k+q) ~ 
-(J4 + J3 + A^i + A^2 + Ar3 + Ar4)] 



1^4 COth(^)(/3 + /4 - Li - L2 - L3 - L4) 



(3.139) 



X4 = |{ (nk+q + rik + 1) [-^^3(/l + /2 + i^l + i^2 + + K^) 

- coth(^)(Ji + J2 - Ml - M2 - M3 - M4) 
-(rik - rik+q) [-^4:{h + /4 + -^vi + L2 + L3 + L4) 

- C0th(^)(J4 + J3 - iVl - iV2 - iV3 - iV4) 

By using the expressions of Appendix E, we obtain the elements in terms of Z**. 
The approximation for the function n(q^, t', t") is given by : 



(3.140) 



or 



^(q^ (',«"> 



(i^A; c/^pi d?p 



(3.141) 



(27r)3 (27r)3 (27r)3 



/**(p + q, -p, k, -q - k, t', to)S**(p, -p - q, Pi, q - Pi, to, t") 

(3.142) 

It remains now to solve the dynamical equations for /** taking into account of the 
loop term J ( |3.118| ). We will do it in the static case. Equations ( p.ll4| )-( |3.117| ) form a 
nonhomogeneous first order differential system. The solutions /**(k, — q — k, pj, — pj + 
q, t',t") are given in Appendix F. The approximation to first order for n(q^,t',t") is 
composed of two terms which are for t' < t" : 



^(q^^^t" 

tr 



i f d^k d^pi d^p 



■X 



2 J (27r)3 (27r)3 (27r)3 

/**(p + q, -p, k, -q - k, t', to)(2o(p) - r)/**(p, -p - q, Pi, q - Pi, t", to)(2o(p + q) + r) 
tr[/**(p + q, -p, k, -q - k, t', t,)X{t,, t")] 

(3.143) 
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The first term corresponds to the approximation one obtains without optimization with 
respect to the initial state (by replacing So by Hq) . When one optimizes with respect to 
initial state, one has to take into account of the second term. We insert in the previous 
formula the expressions obtained to first order for the solutions in the static case. We 
obtain the following result at the first order, for At = t' — t" < : 

n(^)(q2,t',t") = / - — r^^-^ X 

(nk+q + ^k + 1) |^(n-pi+q + rip, + 1) ^2^^Q2 " (^-Pl+Q ^ ^Pl) ^2^^^2 

X ^coth(^^) cos(03At) + i sin(n3At)^ 

-K+q - n^) (^{n-p,+^ + rip, + " ("-Pi+q " ^P^^?]f~^) 

X ^coth(^^) cos{n4At) + i sin(^]4A^)^ 

(3.144) 

with 

^3 = ^^k+q + tUk , fi4 = t^k+q-^k (3.145) 
^7 = ^^-pi+q + t^pi , ^8 = ^^-pi+q - ^^-pi (3.146) 

This is the essential result of this section. One has to add to this formula the term 
( |3.125| ) obtained to lowest order. In the static case, we obtain an approximation for the 
correlation function n(q^,t',t") which depends only on the time difference At = t' — t" . 
If we do not optimize with respect to the initial state, which corresponds to consider only 
the first term in ( 3.143| ), the approximation obtained for n(q^,t',t") contains a spurious 
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dependence on the initial time to which writes : 

p(i)= ' ( ^'^1 ^ X 

2 7 (27r)3 (27r)3 4u;kt^k+qt^pit^-pi+q 

+ -K+q + ^k + l)(n_p^+q + np^ + 1) ^^ _ ^2 m7Coth(^) - i23COth(-^ 

3 7 \ 

X [cosr27(to — t") cosVL'iito — t') + cosVL-jito — t') cosVLziio — t")] 
+ ^(^^k+q + ^k + l)(n_p,+q + np, + 1) ^2 ]_ ^2 (^3 coth(^) - f]7COth(^ 

3 7 \ 

X [sin VtyitQ - t") sin fig (to - t') + sin VL^{tQ - t') sin fig (to - t")] 
-^(^k+q - r2k)(^-pi+q + rip, + 1) ^2 coth(^) - fi4 coth(^) ^ 

X [cosfi7(to — i") cosfi4(to — i') + cosfi7(to — i') cosfi4(to — t")] 
-2K+q-r2k)(^-Pi+q + r^Pi + 1) ^2 _ ^2 I "^7 coth(^) + fi4 coth( — 
X [sin fi7(to - t") sin fi4(to - t') + sin fi7(to - t') sin fi4(to - t")] 
-^(^^k+q + ^^k + l)(n_p,+q - np J ^ ^ (fig coth(^) - fig coth(^ 

3 8 \ 

X [cosfi8(to — t") COS fig (to — t') + COS fig (to — t') COS fig (to — t")] 

-^(^k+q + ^k + l)(r2-pi+q - ^Pi) ^2 _ q2 I ^^3 COth(-— ) -fi8COth( — 

3 8 \ 

X [sin fig (to - t") sin fig (to - t') + sin fig (to - t') sin fig (to - t")] 

+ ^(^k+q - ^k)(r^-pi+q - ^Pl) ^2 1 ^2 (^^8COth(^) - fi4 COth(^)^ 

X [cosfig(to — t") COsfi4(to — t') + COsfig(to — t') COsfi4(to — t")] 

+ ^(^k+q - rik)(n_p,+q - ^pi) ^2 1 ^2 (^fi4COth(^) - fig COth(^)^ 

X [sin fig (to - t") sin fi4(to - t') + sin fig (to - t') sin fi4(to - t")] 

(3.147) 

When we optimize with respect to initial state, the second term of ( |3.143| ) tr[Z**(t', to)X(to, t")] 
cancels exactly the term P^^\ One checks also that at t' = t" = to, this second term is 
vanishing. 

In order to obtain the expression of n(q^, t',t"), we have used a symmetry in the 
exchange of k and — pi of the function to integrate in n(q^, t', t"). It is important to note 
that the function / -^^Sf *(k, — q — k, pi, q — pi, t', t") remains dependent on to. 

For t' > t", expression ( |3.143| ) becomes : 
ii(q,t,t j-^y (2^)3(2,r)3(27r)3'' 

tr [/**(p + q, -p, k, -q - k, t", to)(E:o(p) - r)/**(p, -p - q, pi, q - pi, t', to)(So(p + q) 4 
- tr[/**(p + q, -p, k, -q - k, t", to)X(to, t')] 

('3.148^ 



T 



(3.148) 
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and one shows that one obtains n(^''(q^, At) from ( |3.144| ) by replacing At by —At. 

Similarly to the case of the Green function with two field operators $(x)$(y), we can 
check the fluctuation-dissipation theorem for the Green function with four field operators 
$^(x)$^(y). The expression of n(q^, r) for At > is indeed equivalent to expression 
( gA44l ) of n(q2, -r - for At < 0. 



Let us study the retarded Green function with two field operators x<i'*,<i>*(''^') 
defined by ( |2.120|) when we optimize with respect to the initial state. For t' > t" and with 
the time-ordered product : 



Sf*(k, -q-k,pi,q-pi,t',t") 

For t' > t" and the anti-T product : 
Eff(k, -q-k,pi,q-pi,t',t") 



S**(p, -p - q, Pi, q - Pi, to, t'] 



(2^ 



tr 



/**(p - q, -p, k, -q - k, t", to) 



(3.149) 



(2^ 



tr 



^ (p-q, -P,k, -q-k,t',to) 



S (p, -p-q,pi,q-pi,to,i"' 



(3.150) 

We use formula (|3.134| ) obtained at the first order in perturbation. The retarded Green 
function X'i><i>,'i>'S>(t' ,t") is given by expression ( ^.121| ) obtained by optimizing only the 
dynamics and a term equal to : 



(2^) 
d^p 



+ J ^tr[/**(p-q,-p,k,-q-k,t",to)X(to,t')] 



(3.151) 



where tr[Z**(t', to)X(to, t")] is equal to the term P^^) (eq. (I^M^)). We see from the 
expression of P(^) by interchanging t' and t" that the difference tr[/**(t", to)X(to, t')] - 
tr[/**(t',to)X(to,t")] is vanishing. 

To conclude, for the the linear response formula which involves the retarded Green 
functions, we obtain the same result with or without the optimization with respect to 
initial state. It is sufficient to solve the backward dynamical equations ( p.ll4| - p!ll7| ) in 
real time for Z**. 
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The first order contribution to the polarization function writes : 



n 



(1) 

R 



(rip+q + rip + 1) sin fii At j 



(i^k 1 



[(^q-k + + 1) 



(^q-k - '^k, 



[n 



p+q 



rip) smVL2/S.t 



(27r)3 cuka;q_k^"^^ ' '^'^ ' 

[(?^q-k + ^k + 1)7^2 — 7^ ~ (^q-k - ^'^J 



^4 J 



(27r)3 Wk^^q-k 



with 



^1 — ^i^p+q + ^ip ^2 



P+q 



Coir 



tUq-k + Wk — tUq-k — l-i-'k 



(3.152) 

(3.153) 
(3.154) 



and At = t' - t" . 



By taking its Fourrier transform in time and adding the zero order term and the 
first order term, we obtain : 



n^(q2,cu) = j 



[n 



P+q 
h 



n 



1 

(27r)3 t^ptUp+q 



X 



1 + 



p ^ 
d^k 



in 



4 J (27r)3 CUkCUq-k 
^^2 



+ irjy — Ql 

— [(^q-k + "-k + 1/ ^2 



P+q 

' b 
i + - 



n 



(CU + Z77)2 - fi2 



d^k 



-[(r^q^k + r^k + l)^2_^2 



4 J (27r)3 C^k^^q-k 

(3.155) 

Although the coupling constant appears linearly in the numerator, this formula includes 
nonperturbative contributions since the frequencies involves the self-consistent renormal- 
ized mass m^(/5) eq. ( ^.6(j| ). 

The divergent part of the integrals at high momentum is entirely determined by the 
contribution at T = 0. This one writes : 



n«(q2,^) = I 



d^p 



to, 



P+q 



b 

1 + T 



(27r)3 iOpUJp+d (a; + ir/)^ - (a;p+q + ujp)"^ 

d^k 1 CJq+k + tuu 



X 



(3.156) 



4 J (27r)3 t^k^^q+k (t^p+c 



UJr. 



(We remark that there is no term it] in the denominator of the second integral.) 

At zero temperature, IlR{q^,uj) depends only on the variable 3 = 00"^ — q^. At the 
lowest order, 



n2^(.) 



d^p 
(2^ 



LU- 



P+q 



p^p+q 



{uj + iriY - (cjp+q + ujpY 



(3.157) 
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which is related to the exphcitly Lorentz invariant form : 

r d'^v 1 1 
Uf(s) = i / (3.158) 

where q = {uj,q). We have Rell^^\s) = ReHpis) and Imll^l^\s) = sign{uL}) I mil f{s). rfi is 
the renormahzed self-consistent mass at zero temperature given by eq. ( p.64|) . By using 
for instance the results of Kerman and Lin iTBI , 



where A is a momemtum cut-off, /i is an arbitrary mass scale and the function f{s) is 
defined by : 



f{s) = 2 + \ log ^ ^ = si s>Am^ (3.160) 



/(.) = 2-2W^tan-M-^ si < s < Am' (3.161) 



, /s — 4m^, y/s — 4m^ — -v/i 

= 2 + J log , V 3.162 

V s Vs + 4m2 + 

Adding the zero order term and the first order term, we have at zero temperature : 

nHW^nS'w(i-^.og(^)) (3.163) 

and U^r\s) is given by expression (|3.159| ). We see that at the order we have considered, 
the logarithmic divergence appears squared. By resumming the perturbative serie, we 
obtain : 

We can define a renormalized coupling constant gn^fi) by 

2gi,ili)nf{s)^bna{s) (3.165) 

This renormalization is identical to the one obtained by Kerman and Lin in the symmetric 
phase |jl6 . 
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4 Conclusion 



In this paper, we have shown how to derive variational approximations for two-time cor- 
relation functions in theory. We have studied in particular the two-time correlation 
functions with two field operators $(x) and with four field operators $(x) in the symmetric 
phase. By using the variational principle introduced by Balian and Veneroni ||^, we have 
shown how to take into account correlations both in the dynamics and in the initial state 
while keeping uncorrelated ansatze for the variational objects. The calculation sheme 
we have presented is valid for arbitrary time-dependent problems. As an illustration, 
we have calculated explicitly the two-time correlation functions with two field operators 
and with four field operators at equilibrium in the symmetric phase. The approximation 
we have obtained for the two field operators correlation function is identical in this case 
with the usual mean-field result. This won't be true in the asymmetric phase or for a 
time-dependent problem. The approximation for the two-time and four field operators 
correlation function contains already in the simple case we have considered corrections to 
the standard mean- field result. For the non- retarded correlation functions, it is important 
to optimize with respect to the initial state in order not to obtain a spurious dependence 
on the initial time. However, for the retarded correlation functions, the same result is 
obtained whether or not we optimize with respect to the initial state. Our variational 
approximations for the two-time correlation functions satisfy the fluctuation-dissipation 
theorem. For the four field operators correlation function, we were obliged to solve the 
dynamical equations perturbatively. The zero order approximation was already known 
in the litterature but we have shown here how to derive it variationally. We have cal- 
culated also the first order approximation. At zero temperature, it is possible to resum 
the perturbative serie and to define in this way a renormalized coupling constant. We 
stress that our approximations remain nonperturbative since they involve a self-consistent 
renormalized mass. 

There are several important questions for future works. First, how can one obtain 
variational approximations which include damping. It would be very interesting to ob- 
tain a variational nonperturbative approximation for the viscosity defined in eq. ( 3.86|) . 



Damping effects are probably related to a non symmetrical form in x and y for the matrix 
T defined in eq. ( |2.24| ) to characterize our gaussian variational ansatze and to a choice 
for the operator H different from the Hamiltonian H. A second interesting problem is the 
evolution of a scalar field in self-interaction in a time- dependent metric. The difficulty 
in this case is to solve the backward dynamical equations in order to obtain variational 
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approximations for the two-time correlation functions. 
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Appendix A : Characterization of a Gaussian state 

We call a Gaussian operator an operator which is an exponential of quadratic and 
linear forms of the field operators $(x) and n(x). A Gaussian operator T>{t) is completely 
characterized by the vector a(x, t) and the matrix S(x, y,/!:) defined by the equations 
(^l8D-(g:25l) of the section 2.3. 

Instead of the quantities G, T and 5*, the authors of reference [P] have introduced 
the quantities G^, and ^, ^ being the degree of mixing (for a pur state, ,^ = 1) : 

Tr(P(t)$(x)<|.(y)) = (Gf (1 - O^^Gf ) (x, y, t) , (A.l) 

rr(P(t)n(x)n(y)) = 1 {Gf\l + OG^'/') (x, y, t)+4 (s^Cf (1 - O^'^f S^) (x, y, t) 

(A.2) 

Tr(P(t)$(x)n(y)) = 6%^ - y) + 2 {Gf{l - O^'^f S^) (x, y, t) . (A.3) 
It is usefuU to relate our notations to those of the authors of 0. For a pure state 

(e = o), 

T(x, y , t) = 2 (GS + SG) (x, y , t) , (A.4) 

5(x, y, t) = ^ G-' (x, y, t) + 4 (SGS) (x, y, t) . (A.5) 
For an uniform configuration, we have the following relations : 

G{p)=G^{p){l-ap)r\ (A.6) 

^(P) = \ Gfip) (1 + e(p)) + 4S|(P)G^(P) ' 

T(p)=4G5(p)^-i^S5(p) , (A.8) 
and the Heisenberg invariant is related to the degree of mixing according to : 
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Appendix B : Expression of the von Neuman entropy for bosons in the 
Hartree-Bogoliubov approximation 

The von-Neuman entropy is given as a function of the density operator T) by : 



TrT) log T) 

S = +\ogTrV . 

TrV ^ 



(B.l) 



In the HB approximation, V is the exponential of linear and quadratic forms of the field 
operator $ and 11. The entropy density then writes in terms of the matrix H : 



S = I 



(2^)^ 4 



^ trK-il-T^-Eil-r)-!] lo ^^^^ 



(B.2) 



where r is the 2x2 matrix 



r 

-1 



It can also be written in terms of the Heisenberg invariant I ( p. 261) : 



S 



(27r)= 



(B.3) 



(B.4) 



For the static HB solution characterized by S'^, we have : / = 
easy to convince oneself that expression ( |B.4| ) is in this case identical to 



/dp 
j^-^ [{rip + 1) log(np + 1) - rip log rip] 



^0 T| j 



(B.5) 



where the occupation number Up is related to the degree of mixing ^(p) according to : 



2rip + 1 



i + e(p) 

\ 1 - e(p) 



Hp) ■ 



(B.6) 



By using the parametrization ^(p) = ^^^^^^ , we have also : 2rip + 1 = coth(^^). In the 



^'^ theory, Up = \J-goip) and 

^o(p) = 
The static solution satisfies : 



P'^ + ml + ^^l + ^Go{x, x) 



(B.7) 



2?i(p) ^7o(p)=H°2(p) , 



(B.8) 
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and 



s;,(p) s^,(p) = 



i + e(p) 
i-e(p) 



In terms of the frequency ujp : 



S?i(p) = COth(^(^p) 
-22{P) = -^P COth(^Wp) 



/(p) = coth^(^tc;. 



(B.9) 

(B.IO) 

(B.ll) 
(B.12) 



The matrix of the second derivatives of the entropy evaluated for the static solution 
So is given by the following elements : 

,2 



asii(p)(9Sii(p) 



asii(p)as22(p 



UJZ. 



1 



X 



+ 



coth X coth^ X — 1 

1 1 



X 



4 L coth a; coth^x — 1 



5522 (P) 5^22 (p) 



X 



1 



95i2(p)9Si2(p) 



coth a; coth^ x — 1 

X 



cothx 



(B.13) 

(B.14) 
(B.15) 
(B.16) 



where we have introduced the notation x — j3ujp/2. 



Appendix C : Parametrization of the product of two Gaussians 

In this appendix, we give the expressions of the expectation values a^, S;, and Sc 
in terms of a^, «„, and 'E^- 



~r ~r ^d 



where r is 2 x 2 matrix 



1 



r 



.-1 oy 

For the matrices S, we have the following relations : 

Sfe - T = (S„ - r) ^ (Sd - r) , 



(C.l) 
(C.2) 

(C.3) 



(C.4) 
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r = (S, - r) (S„ - r) . (C.5) 



Lb + T = (S<i + r) ^ ^ (5a + r) , (C.6) 
:e + T = (5„ + T) (5, + t) . (C.7) 



Appendix D 

In this appendix, we give the evolution equations for na, aa, Ua, aa, Sq 



(D.5) 



(D.2) 



2i = - [{Ed + r) {He - I'k) (Sd - r) - (S^ - r) Tib (5d + r)] r a^-e 
+Sd (wc - Wb) + r {wc + Wft) - (Sd + r) 

i = - [(S, + r) (Tie - X^) (S, - r) - (S, - r) H, (S^ + r)] (D.3) 



2i d„ = [(S„ - r) (K - ^f'^) (5a + r) - (S„ + r) Hb (S„ - r)] r a^-b 
+Sa (wb - Wc) + T (wc + Wb) - (S„ - r) ty^ 

i 5„ = [{Ea - r) (Tic - Tk) (2a + r) - (S, + r) Tib (S^ - r)] (D.6) 



The expressions for the vector w and the matrices Ti and I are the following : 
W6(x, t)i = ^ f ^J' = -/^(x, t) , Wb{x, t)2 = i 'l^ = ^7r6(x, , (D.7) 

Tib(x, y, t),, = -2 ^ ^ ^ , (D.8) 

^6(x,y,i)ii = -^^^^-^ = ^i?6(x,y,i) , (D.9) 
dGfe(y,x,t) 2 
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(D.12) 



(D.13) 



The quantities F which appear in equations (D.l) and (D.3) are defined by : 

If 

+- J (fx(fy tr ['^^(x, y, t){Eh - + 2ab-dab-d){y, x, t)] 

FrI = / (A^c(x, t) - WKci'^, t)ac-d{y^, t)) 
I r 

+- J d^xd^y tr [Xkc(x, y, t){Ec - + 2ac-dac-d){y, x, t)] 
For the theory, we have : 

/6(x, t) = - ^-A + ml + ^ ifili^, t) + ^ Gb{x, X, t) j <^,(x, t) , (D.14) 

^6(x, y,t)^- (^-A + ml + ^ t) + ^ ^^(x, x, t)^ 5(x - y) . (D.15) 
From the variations of the source term K^, we obtain : 

<(x,t)i = J*(x,t)+2 y^^xa (j**(x,X2,t)(^c(x2,t) + J*"(x,X2,t)7re(x2,t)) , (D.18) 

«;^(X, i)2 = i J"(X, + 2 i y d^X2 (J™(x, X2, (/?c(x2, t) + J™(x, X2, 7rc(x2, t)) , 

(D.19) 

T^(x,y,i)n = -J**(x,y,i) , (D.20) 

J^(x, y, t)i2 = -2i J*"(x, y, t) , (D.21) 

X^(x,y,t)22 = J™(x,y,t) . (D.22) 

The expressions for the matrices t,T,r and R which appear in the dynamical equa- 
tions for the two-time correlation functions are the following : 

tn(x,y,i) = -^^°)(x,y,i) , t22(x, y, = -(^(x - y) , (D.23) 

ri,n(x,y,z,t) = -^(^(°)(x,i)5(x-y)53(x-z) , T^jk = , (D.24) 

rn,i(x,y,z,t) =ri,n(x,y,z,i) (D.25) 

i?ii,n(x, y, z, u,t) = ^ 6\k - y) ^^(x - z) 5(x - u) . (D.26) 
The other elements are equal to zero. 
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We use the following definitions for the Fourrier transform : 



(2^'"(2^ exp(i^Pj- .Xj) iy(pi,...,p„) , (D.27) 
Vr(pi, p„) = j d^xi...(fxn exp(-i Yl Pi • ^i) x„) , (D.28) 



\3 x3/ 



iy(pi, p2, p„) = {2tiY 5^(pi + p2 + ... + p„) iy(pi, P2, Pn) ■ (D.29) 
So for the two-point function /3* : 

/?f(x,xi,t,t") = j ^e^P-(''i-)/?f(p,t,t") . (D.30) 



Appendix E 

We give in this appendix the expressions of the integrals In, Jn^ Kn, Ln, M^, Nn, 
n = 1,2,3,4 in termes of the constants a,b,c,d introduced in ( p.l29|) and of We 



introduce the frequencies 



We define : 



_ f dH a ^ r dH c 

r d^l a r d^l c 



^3/ 


a 




(27r)3 






dH 


a 




(27r)3 


ni- 








5 


(27r)3 


m- 




rf3/ 


bQ 


5 


(27r)3 


ni- 





d^l c , , /■ (^3/ (i 



/■ d I d 

(2vr)3 - nl ' = y (2vr)3fii-ni^^"^^^^^) ^^'^^ 

The integrals L„ are obtained from the integrals Kn by replacing ^3 by ^4. 
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The integrals Nn are obtained from the integrals M„ by replacing Q3 by Q4. 



The expressions of sums of integrals in terms of If,^ are usefuU : 



(27r)3 



ni + 1 



^3 ^5 



1 



1 



-/ff +z2*'')coth(^) + -/r 



1 



^3 ~ ^6 



12 ^ 21 

(^/f * - itr) coth(^) + i/f * - J- 

a;ia;i+q z a;i ^i+q 



"''21 



(E.ll) 

/s + I A is obtained from Ji + /2 by replacing by O4. The arguments of the functions 
*■* in these integrals are If^^ 



If^ in these integrals are — q — 1, pi, q — pi, i", to). 



Jl + J2= V 



(2vr)^ 



■X 



^2_^2 (^1+q + ^l + l) 



^ /fi* + /|*) + (-^f2* + ^/2T)coth(^; 



■i '4><I> 74><I>\ I /_|_7 

11 '22 J "T I . '12 



^1^1+q ^1 

J3 + J4 is obtained from Ji + J2 by replacing Q3 by Q4. 



^^l+q 

^21 ) COth(— -) 



2) 



dH 



{2nf 
ni+q + ni + 1 

^3 ~ ^5 



X 



-/ff + /f*)coth(^)-(-/f* 



1 



1 



7<1><I>\ 

"'■21 J 



'^1+q 



;^/f * - itf) coth(^) + l/f^* .21 

wiWi+q z a;i '^l+q 

(E.13) 



Li + L2 + L3 + L4 is obtained from Ki-\- K2 + Kz + by replacing VL^ by O4. 



Ml + M2 + M3 + M4 ^ -V j 



■X 



^3 ~ ^5 



dH 

(ni+q + ni + 1) 



1 



22 J 



199 , 7 
I'll 



1 /$* I 
'19 "T 



6i;i 



— ^(m+q-m) 



<^l+q 



-i2i ) coth(- 



6i;ia;i+q 

'^—ltl-lt2) + {-lt2-^ltl)^Oih{ ^ 



2 



(E.14) 



Ni + N2 + N2, + N4 is obtained from Mi + M2 + M3 + M4 by replacing 1^3 by ^l^. 



Appendix F 



57 



In this appendix, we give the solutions of the backward dynamical equations for I 



ij 



at the first order . By noting : 

^3 = (^k+q + <^k , ^^4 <^k+q " <^k 



^5 — '^Pi-q + ^Pi ) ^6 — CO 



Pi-q ^pi 



4a;p.a;q_p, 



(F.l) 
(F.2) 
(F.3) 



^n*(k, -q - k, p^, -p^ + q, t' , t)) = 

-[S^ik + Pi) + 5^{-q - k + Pi)] cosa;q+k(i - t') cos ujk{t - t') 



2' 
F 

-^K-Pi + ^Pi + 1) 



cos 1^3 (t - t') + 



4 — ^2,5 



cos 04(t — t') 



+ ^(^q-Pi - ^Pi 



^,^C0sQ3(t-0 + ^^ 



cos Q4{t — t') 



^5 +_i^]cosn5(^-0 



■-(nq_p. -np.) 



6 , ^6 



+ 



i^Q ^^/l ^^fl 



COS VL^it — t') 



(F.4) 



^12 (k, -q - k, Pi, -p, + q, t', t)) = 

[5^{k + Pi) + S^{-q - k + Pi)] sina;q+k(i - t') cosu}k{t - t') 



2ia;k+q 
F 

2ia;k+q 
F 

4 



n, 



q-Pi 



+ «Pi + 1) 



^5 / /\ ^5 

^^^sinl]3(t-t) + ^^^ 



sin Q^it-t') 



2ia;k+q 



sinQ3(i-i') + 



sin ^^{t-t') 



+ 



Q4 



i^q 



sin f]5(t-t') 



^fi*(k, -q - k. Pi, -Pi + q, t', t)) = 



sin (t — t') 



-^[^'■^(k + Pi) + S^{-q - k + Pi)] cosWq+k(t - t') smu!k{t - t' 
F 

J— (nq_p. +np. + 1) 



2ia;k 

+77— K-p, -np.) 



sin n4{t-t') 



F , 

+77— K-p, +np. + 1) 



sinl]3(t-t') - 



Q4 



sin f]4(t-t') 



2ia;k 
2ia;k 






^6 ^4 ^6 



sinQ5(^-0 
I sinOel^-^') 



(F.5) 



(F.6) 
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^22 (k, -q - k, Pi, -Pi + q, t', t)) = 

+- [6%k + p,) + (53(-q - k + p,)] sin Wq+k(t - t') sinco'k(t - t') 

^ COS Qs{t - t') - „2 ^r>2 ^4:{t - t') 



2a;ka;k+q 
F 

+ 

+ ^ in, 
2a;k<^k+c 



nq-p, + Up. + 1) 



^3 ^5 



q-Pi 



Pi J 



cosflsit-t')- ^2_o2 



COS 04(i — t'] 



(F.7) 



q-Pi 



Sl^kf-^k+q 



,^q-pi 



+ np, + 11 



O5 






cos fl5{t — t') 



^^Pl ^^f, 



cos Q6(^ ~ t') 
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